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ABSTRACT 


An  algorithm  pertinent  to  pattern  classification  and 
cluster  analysis  of  n-dimensi onal  (n>2)  data  is  presented. 
The  algorithm  simplifies  the  dimensionality  of  the  sample  set 
by  a  perpendicular  projection  of  each  point  onto  a  plane  which 
is  determined  by  two  orthonormal  vectors.   The  operator  views 
the  plane  on  a  computer  graphics  screen,  selects  a  point  on 
the  screen  with  a  lightpen,  and  attempts  to  achieve  separation 
of  classes  by  choosing  a  "direction  of  movement."   Each  itera- 
tion, a  movement  of  a  single  point  in  a  certain  direction,  is 
actually  a  reorientation  of  the  plane  into  which  the  data 
points  have  been  projected.   These  iterations  continue  until 
the  operator  decides  successful  separation  has  been  accomplished 
The  method  is  neither  restricted  by  statistical  distribution 
assumptions   nor  the  dimensions  of  the  sample  set.   From  the 
investigation  of  seven  case  studies,  this  algorithm  emerges 
as  a  tremendous  assistance  in  the  analysis  of  multidimensional 
data. 
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I.   INTRODUCTION 


The  algorithm  presented  in  this  paper  is  an  improvement 
on  the  separation  algorithm  presented  in  the  thesis  by  Gilbert 
Lauzon  [1].   He  developed  a  separation  algorithm  to  use  in 
pattern  classification  and  cluster  analysis  of  n-dimensi onal 
data. 

The  algorithm  was  written  in  Fortran  IV  for  the  XDS  9300 
computer  installation  at  the  Naval  Postgraduate  School.   The 
Adage  Graphics  and  Text  (AGT)  system  is  linked  directly  to 
the  digital  computer.   The  Fortran  callable  graphics  sub- 
routines in  the  XDS  9300  allow  for  communications  between  the 
two  systems.   Also,  by  using  the  library  program  "MAD"  of  the 
AGT  system,  the  lightpen  and  sixteen  function  switches  permit 
an  alteration  of  the  display  on  the  graphics  screen. 

The  computer  program  projects  each  point  into  a  plane 
which  decreases  the  dimensionality  of  the  analysis  from  the 
original   n  (>2)  to  2.   The  points  on  the  plane  are  presnted 
to  the  operator  on  a  cathode  ray  tube.   The  user  reorients 
the  plane  on  each  iteration  by  selecting  a  point  and  choosing 
a  "direction  of  movement."   Using  the  XDS  9300/AGT-10  system, 
the  operator  makes  the  decision  whether  to  proceed  with  a 
different  orientation  of  the  plane  or  terminate  the  analysis. 
The  human  decision  is  based  on  the  pattern  of  points  presented 
on  the  screen.   The  interaction  of  the  man-machine  system 
is  emphasized. 


Many  classification  problems  can  be  approached  in  a 
graphical  diagram.   A  two-dimensional  plot  of  the  features  is 
drawn.   From  this,  the  human  must  suggest  a  rule  that  partitions 
the  feature  space  into  regions,  one  for  each  category.   The 
rule  can  be  a  series  of  straight  or  curved  lines  depending  on 
the  number  of  categories.   Since  the  decision  is  based  on  the 
perception  of  the  operator,  the  rule  is  quite  arbitrary. 
Therefore,  the  analyst  must  insure  on  a  proper  number  of 
samples   in  each  category.   In  addition,  the  human  tries  to 
determine  the  rule  which  according  to  his  subjective  judgement 
minimizes  the  probability  of  error.   In  essence,  this  means 
that  the  probability  of  a  random  observation  from  a  specific 
category  falling  outside  its  respective  region  defined  by  the 
rule  should  be  made  as  small  as  possible.   The  preceding 
analogy  is  the  basic  characteristics  that  prevail  in  the 
separation  algorithm.   The  graphical  diagram  is  the  display 
screen  with  the  projections  from  the  computations  of  the  digital 
computer.   The  rule  is  solely  determined  by  the  operator  from 
the  video  presentation  of  the  graphical  diagram. 

The  first  modification  incorporated  into  G.  Lauzon's 
algorithm  was  to  insure  that  the  A  and  B  vectors  were  ortho- 
normal  on  each  iteration  through  the  computer  program. 
Basically,  this  means  that  the  two  vectors  are  perpendicular 
and  each  vector  has  a  length  of  one. 

Secondly,  the  existing  "X-Y"  video  display  was  changed  to 
a  more  discernible  pair  of  symobls  consisting  of  the  "X-0" 


scheme.   This  improvement  greatly  enhanced  the  ability  of  the 
operator  to  visually  identify  the  two  groups  on  the  graphics 
screen . 

Thirdly,  supplementary  video  information  was  added  to  the 
display  on  the  screen.   Each  iteration  of  the  algorithm 
involves  the  selection  of  a  point  and  a  "direction  of  movement." 
These  operations  are  controlled  by  a  sequence  of  steps 
utilizing  the  sixteen  function  switches  and  the  lightpen.   The 
function  switches  have  been  numbered  left  to  right  and  top  to 
bottom  for  easy  reference.   A  record  of  the  use  of  these 
function  switches  is  most  beneficial  to  the  operator.   So, 
the  numbers  of  two  function  switches  have  been  displayed  in  the 
upper  left  hand  corner  of  the  screen.   The  one  number  records 
the  function  switch  corresponding  to  the  type  of  symbol 
selected  and  the  second  keeps  a  record  of  the  other  function 
switches.   The  modification  was  incorporated  to  merely  assist 
the  operator  and  did  not  contribute  to  improving  the  separation 
algorithm. 


II.   DISCUSSION  OF  THE  PREVIOUS  ALGORITHM 

The  previous  separation  algorithm  was  presented  in  the 
thesis  by  Lt.  Gilbert  Lauzon  [1].   The  algorithm  is  applicable 
to  use  in  pattern  classification  and  cluster  analysis  of 
n-dimensional  data.   The  computer  graphics  system  projects 
each  point  of  a  data  set  into  a  plane.   The  projection  is 
not  necessarily  orthonormal .   The  method  reduces  the  dimen- 
sionality of  the  problem  from  n  to  2. 

First,  a  description  of  the  algorithm  will  be  presented 
for  a  pattern  classification  problem.   Assume  there  are  two 
classes  of  objects,  P-,  and  P2-   These  two  groups  are  comprised 
of  L-dimensional  vectors  of  numbers.   These  numbers  may  be 
measurements  or  merely  some  assigned  values  of  nature.   The 
goal  of  the  analysis  is  to  determine  the  existence  of  a 
partition  of  the  two  groups.   Once  a  separation  has  been 
achieved  and  a  discriminant  function  has  been  determined,  it 
can  be  used  to  classify  new  objects  by  the  rule: 

f(Z)  >  0  -*■  Z  is  a  member  of  the  set  P,  ; 
f(Z)  <  0  ->  Z  is  a  member  of  the  set  Pp. 


The  two  groups  are  considered  separable  if  an  appropriate 
discriminant  function  can  be  specified.   Otherwise,  the  classes 
are  inseparable  and  no  discriminant  function  can  be  assigned. 


The  L-dimensional  vectors  of  the  measurements  on  objects 
from  the  classes  P,  and  P2  will  be  denoted  by  X  and  Y, 
respectively.   The  algorithm  commences  by  calculating  the 
average  value,  DAV ,  of  each  of  the  L-components  of  the  data 
set. 


DAV(K)  =  ( 


1  . 


M 


N 


M  +  N 


)  (  Z  Y(I,K)  +   Z  X(J,K))  ;  K  =  1  ,...  ,L 


1  =  1 


J  =  l 


Here  M  and  N  are  the  number  of  vectors  in  P,  and  P?  ,  respec- 
tively.  Using  the  L-component  vector  DAV,  the  data  points  are 
centered  as  follows: 


X'(J,K)  =  X(J,K)  -  DAV(K)  ;  J=1,...,M  and  K=1,...,L 
Y'(I,K)  =  Y(I,K)  -  DAV(K)  ;  1  =  1 N  and  K=1,...,L 

Two  arbitrary  non-zero  vectors,  A  and  B,  in  Euclidean  L-space 
are  chosen.   A  linear  transformation i s  computed  for  each  point 
using  the  two  vectors.   This  amounts  to  projecting  each  vector 
to  a  point  in  Euclidean  2-space  with  coordinates 


TAjB  (Z)  =  (A-Z,  B-Z) 

L  L 

where  A»Z  =   Z  avzVi  B*Z  =   Z  bvzVi    1    is  an  arbitrary  point 

K=l  K  K        K=l  K  K 

from  the  X  or  Y  set.  The  coordinates  are  scaled  so  the  points 

can  be  displayed  on  the  AGT  screen.   The  procedure  determines 

the  maximum  absolute  value  of  both  coordinates  as  shown: 


TMAX  =  max  (max  | A  -  X |  ,  max  | A  - Y |  ) 

J  I  I  =  1  ,.  ..  ,N 

J  =  1 ,. . .  ,M 
SMAX  =  max  (max  |B-X|  ,  max  | B - Y |  ) 
J  I 

Each  coordinate  is  divided  by  the  respective  maximum  value 

TA  B  (Z)  =  (A-Z/TMAX,  B-Z/SMAX) 

At  this  point,  the  video  display  will  appear  depicting  the 
projections  of  the  X  and  Y  points.  The  operator  selects  a 
point  and  chooses  a  "direction  of  movement."  For  example, 
if  the  tenth  X  point  is  picked  and  a  move  .right  is  chosen, 
a  new  vector,  A',  is  computed. 

A' (K)  =  A(K)  +  X(10,K)  ;  K  =  1,...,L 

The  new  values  of  the  A  and  B  vectors  are  dependent  on  the 
"direction  of  movement"  and  the  selected  point  as  follows: 

Direction  of  Move        New  A        New  B 


RIGHT 

A  +  Z 

B 

LEFT 

A  -  Z 

B 

UP 

A 

B  +  Z 

DOWN 

A 

B  -  Z 

The  computer  returns  to  the  linear  transformation  calculations 
and  repeats  the  procedure  using  the  new  vectors. 

The  goal  of  the  procedure  is  a  separation  of  X's  and  Y's 
into  disjoint  areas  of  the  plane  (screen).   The  operator 
proceeds  using  the  lightpen  and  sixteen  function  switches  until 


an  appropriate  separation  is  attained  or  discontinues  the 
analysis.   If  the  two  classes  have  been  successfully  separated, 
the  discriminant  function  f(Z)  can  be  specified.   For  instance, 
if  the  X's  and  Y's  on  the  screen  can  be  separated  by  a  straight 
line,  the  discriminant  function  would  be  of  the  form: 

f(Z)  =  d  (A-Z/TMAX)  +  e  (B-Z/SMAX)  +  f 

where  d,  e,  and  f  are  scalars.   The  function  need  not  be 
unique.   The  pattern  classification  is  completed  once  the 
discriminant  function  has  been  determined. 

The  separation  algorithm  for  the  cluster  analysis  problem 
is  yery   much  like  that  in  the  pattern  classification.   Normally, 
the  set  of  data  is  not  broken  down  into  two  classes.   The 
entire  set  is  presented  as  a  set  of  Y's.   The  concept  of  this 
method  is  the  determination  of  groupings  from  the  sample  set. 
It  is  normally  used  on  data  where  the  goal  of  the  analysis  is 
to  detect  clusters.   The  final  groupings  of  the  data  points 
into  the  appropriate  number  of  disjoint  subsets  is  the  decision 
of  the  operator.   Each  of  the  subsets  represents  a  cluster. 
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III.   DISCUSSION  OF  THE  NEW  ALGORITHM 

As  mentioned  in  Section  I,  the  separation  algorithm  was 
modified  using  orthonormal  vectors,  A  and  B.   This  modifica- 
tion is  important  for  both  the  pattern  classification  and  the 
cluster  analysis  problems.   The  following  paragraphs  will 
describe  in  detail'  the  concept  of  the  orthonormal  projection 
al gori  thm. 

First,  let  us  assume  the  analysis  is  conducted  on  a  pattern 
classification  data  set.   The  data  has  been  divided  into  N 
vectors  in  P-,  and  M  vectors  in  P2-   The  P-,  set  is  read  by  the 
computer  as  the  Y  points  and  the  P~  set  as  the  X  points.   The 
data  points  have  been  centered  by  calculating  the  L-dimensi oned 
vector,  DAV ,  by  averaging  each  component  of  the  entire  data 
set. 

Next,  the  operator  must  select  the  A  and  B  vectors  to  be 
inserted  into  the  computer  memory.   If  the  inner  product 
between  two  vectors  is  zero,  then  the  two  are  perpendicular. 
The  length  of  any  vector  is  determined  by  taking  the  square 
root  of  the  inner  product  with  itself.   For  example,  the  length 
of  a  vector  A,  denoted  by   |A| I ,  is  given  by 


||A||  =  (A- AT2  =  (  Z  a/) 


K=l 


K 


Now  the  linear  transformation  is  performed  on  each  point 
as  follows: 


TA  B  (Z)  =  (A-Z,  B-Z);  A-Z  =   Z  a^,  B-Z  =  J-^K^K 
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L 

Z 

K=l 


This  transformation  is  an  orthonormal  projection  of  the  point 
Z  (an  arbitrary  X  or  Y)  into  the  plane  defined  by  the  ortho- 
normal  vectors,  A  and  B.   The  two  sets  of  projections  are 
presented  to  the  operator  by  the  video  display  on  the  screen. 

The  scaling  method  described  in  Section  II  is  performed 
on  the  data  points  in  order  for  them  to  be  displayed  on  the 
screen.   By  using  only  one  scale  factor,  the  same  perspective 
of  the  two  projections  will  be  retained.   This  scaling  term, 
VMAX,  is  determined  as  follows: 

VMAX  =  max  (  max  ( | A-X |  ,  | B-X | ) ,  max  ( |  A- Y |  ,  | B- Y | )  ) 
J  I 

where  J  =  1,...,M  and  I  =  1,...,N.   The  final  values  of  the 

coordinates  from  the  linear  transformation  will  be 


TA  B  (Z)  =  (A-Z/VMAX,  B-Z/VMAX) 


The  first  iteration  of  the  algorithm  commences  when  the 
operator  picks  a  point  on  the  screen  and  chooses  a  "direction 
of  movement."   Actually,  the  computations  do  not  "move"  the 
point  selected  but  reorients  the  plane  upon  which  the  entire 
data  set  is  projected.   The  reorientation  of  plane  is  accom- 
plished by  computing  two  new  orthonormal  vectors.   Suppose 
that  the  user  has  selected  an  arbitrary  point  Z  and  desires 
to  move  it  right.   The  first  step  in  the  computation  is  that 
of  recomputing  the  A  vector, 

A'(K)  =  A(K)  +  Z(K)  ;  K  =  1  ,...  ,L 
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For  the  projection  of  Z  to  shift  horizontally  to  the  right, 
the  new  vector,  B',  should  be  the  same  as  B.   But  this  is 
impossible  if  A'  and  B'  are  to  be  perpendicular.   Hence,  the 
projection  of  Z  shifts  also  up  and  down.   We  want  this  vertical 
shift  to  be  minimal.   It  is  best  to  follow  the  remainder  of 
the  concept  to  finding  such  a  B  vector  with  mathematical 
expressions . 

The  appropriate  symbols  will  be  defined  as: 

A  and  B   =  old  orthonormal  vectors 
A'  and  B'  =  new  orthonormal  vectors 
| A » B  |     =  absolute  value  of  the  inner  product 
of  the  vectors  A  and  B 

The  new  vector  A'  must  be  changed  so  it  remains  constant  with 
a  unit  length.   This  is  done  by  dividing  the  vector  (A  +  Z) 
by  its  length  as  such: 

A'  =  (A  +  Z )/ | | (A  +  Z)|| 

This  is  the  new  vector  which  will  be  used  in  recomputa ti on  of 
the  new  projections. 

The  new  B  vector  will  be  determined  using  Lagrangian 
minimization.   This  technique  is  applicable  because  the  require 
ment  to  determine  the  minimum  value  of  the  objective  function, 
| Z - ( B '  -  B)|,  has  the  two  restrictions  imposed  on  it.   Now, 
to  find  the  B1  which  minimizes  the  difference  between  Z«B'  and 
Z*B,  the  following  problem  will  be  solved: 
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Minimi  ze 
subject  to 


Z-(B'  -  B) 


B' -A'  =  0 
2 


MB 


=  1 


(1) 
(2) 
(3) 


Using  the  Lagrangian  method  to  find  the  solution,  let 


F  =  +  Z-(B'-B)  -  X1  (B'-A')  -  A2  (|  |B'|  |2  -  U 


(4) 


where  X,  and  X?  are  the  Lagrangian  multipliers.   Taking  the 
first  derivative  with  respect  to  B'  and  setting  the  result 
equal  to  zero: 


|jp  =  +  Z  -  XTA'  -  2A2B'  =  0 


(5) 


Using  equations  (2)  and  (5),  the  value  for  A-,  is 


X1  =  +  Z-A' 


(6) 


This  result  is  true  if  and  only  if  X2  has  a  non-zero  value. 
Using  equations  (3)  and  (5),  the  value  for  X2  can  be  determined 


as  : 


X2  "  t 


1    -  (Z-A1)  A' 


(7) 


Substituting  the  values  of  X,  and  X2  into  equation  (5),  the 
final  solution  to  the  minimization  problem  is: 


R'  -  +  ,(Z  -  (Z-A1)  A') 
B   '  -  I  |(Z  -  (Z-A1)  A') 


(8) 


The  final  step  is  the  determination  of  the  sign  to  use  for  the 
B'  vector.   By  using  the  inner  product  of  Z  with  both  sides  of 
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equation  (8)  and  the  second  derivative  of  equation  (4),  the 
following  statements  can  be  shown  to  be  true: 

(1)  Use  +  sign  if  sign  (Z-B')  f    sign  (Z  -  (Z-A')) 

(2)  Use  -  sign  if  sign  (Z-B')  =  sign  (Z  -  (Z-A')) 

The  same  reasoning  can  be  extended  to  "moving"  an  arbitrary 
point  Z  to  the  left.   However,  in  selecting  a  point  Z  and 
"moving"  it  up,  the  new  vector  B'  will  be  determined  first. 
This  time  the  minimization  problem  will  be: 

Minimize        | Z  •  (A  '  -  A) | 
Subject  to      B'  -A'  =  0 

IIA'M2  =  1 

The  following  table  gives  a  guide  to  the  computations 
of  the  new  vectors  using  an  arbitrary  point  Z. 

Direction  of  Move         New  A  New  B 

RIGHT  (A  +  Z)/| |A  +  Z| |  B* 

LEFT  (A  -  Z(/| |A  -  Z| |  B* 

UP  A*  (B  +  Z)/|  |B  +  Z|  | 

DOWN  A*  (B  -  Z)/|  |B  -  Z|  | 

*Denotes  the  use  of  the  minimization  technique. 
These  new  A  and  B  vectors  are  used  in  the  recomputati on  of  the 

linear  transformation  T  on  the  data  set.   After  the  projections 

have  been  calculated,  a  new  display  of  the  plane  depicting  the 

two  subsets  will  appear  on  the  screen.   The  operator  can  make 
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another  selection  or  terminate  the  analysis  depending  on  the 
presentation  of  the  points. 

The  "X-Y"  letters  in  the  old  routine  were  changed  to  the 
nX-0"  symbols.   This  was  done  to  make  the  presentation  easier 
for  the  user  to  discriminate  the  two  data  sets.   The  similarity 
between  an  X  and  Y  made  it  perplexing  when  the  two  different 
symbols  were  in  close  proximity.   The  two  letters  are  sketched 
on  the  screen  through  a  set  of  draw  and  move  instructions. 
These  computer  instructions  are  stored  in  a  vector  called  IPACK 

Video  information  of  the  function  switches  was  also 
incorporated  into  the  computer  program.   Now,  when  the  operator 
depresses  one  of  the  sixteen  function  switches,  the  number  of 
the  switch  is  displayed  in  the  upper  left  corner  of  the  screen. 
The  number  in  the  first  row  indicates  the  function  switch  for 
selecting  either  an  X  or  a  Y  point.   The  value  of  these 
commands  will  be  an  11  or  3,  respectively.   The  remaining 
function  switches  will  be  indicated  on  the  second  row.   This 
additional  video  display  is  controlled  by  a  series  of  ENCODE 
and  CALL  TEXTO  instructions  in  the  computer  program.   The  role 
of  the  function  switches  is  to  communicate  between  the  Fortran 
graphics  package  "MAD"  and  the  9300  computer.   The  display  of 
the  function  switch  numbers  assist  the  operator  in  efficient 
execution  of  the  program. 
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IV.   ANALYSIS  OF  THE  SEVEN  CASES 

The  operation  of  the  separation  algorithm  using  the 
orthonormal  vectors  was  first  examined  using  four  sets  of 
artificially  generated  coordinates.   The  first  set  was  just 
to  test  the  modified  program.   The  remaining  three  were 
generated  for  the  sake  of  testing  the  applicability  of  the 
method  to  pattern  classification.   Each  set  of  points  filled 
various  geometrical  shapes.   This  quality  allowed  for  ensuring 
proper  appraisal  of  the  procedure  using  mul ti -dimensi onal  data 
Finally,  three  cases  of  measurements  from  "real  world"  applica 
tions  were  analyzed  to  evaluate  the  performance  of  the  new 
algorithm.   The  seven  cases  that  will  be  discussed  are: 

TESTS: 

(1)  One  three-dimensional  pyramid 

(2)  Two  three-dimensional  cubes 

(3)  Two  six-dimensional  cubes 

(4)  One  four-dimensional  shell  enclosing  a  four-space 
sphere 

APPLICATIONS: 

(5)  Iris  flowers 

(6)  Fossils 

(7)  Coronary  heart  disease 


17 


Prior  to  exemplifying  on  each  case,  a  few  comments  are 
appropriate.   The  output  from  the  computer  printer  comprises 
a  sixty-one  by  thirty-seven  dimensional  matrix.   The  printout 
is  a  six  by  six  inch  square.   The  outline  of  the  square  is 
shown  by  asterisks  and  the  coordinate  axes  are  indicated  by 
four  zeroes.   The  symbols  denoting  the  data  points  on  the 
printout  correspond  to  the  plot  on  the  AGT  graphics  screen. 
The  points  will  appear  as  an  X,  0,  8,  E,  or  J.   The  "X"  and 
"0"  symbols  match  the  positions  of  the  X  and  Y  projections, 
respectively.  '  An  "8"  symbol  represents  two  or  more  Y's  in  a 
cell;  and  "E"  refers  to  two  or  more  X's.   A  cell  denotes  a 
matrix  element  of  the  printout  array.   The  matrix  is  stored 
in  the  computer  memory.   For  example,  PIC  (31,18)  in  the 
computer  program  would  refer  to  the  symbol  in  the  cell  occupying 
row  31  and  column  18  of  this  matrix.   Lastly,  the  "J"  letter 
indicates  a  mixture  of  X's  and  Y's  occupying  the  same  cell. 
The  coordinate  axes  have  been  drawn  in  by  hand.   The  encircled 
points  signify  the  X  or  Y  selected  in  the  next  iteration  and 
the  arrow  represents  the  direction  of  movement. 

For  each  case,  the  initial  values  of  the  A  and  B  were 
extracted  from  the  table  in  Appendix  IV.   This  table  represents 
only  one  group  of  an  infinite  set  of  possibilities.   This 
particular  group  was  chosen  because  the  absolute  values  in  any 
dimension  are  nearly  equal.   The  author  felt  the  use  of  nearly 
equal  components  would  initially  bring  all  measurements  into 
the  projection  by  the  same  relative  proportion.   It  was  thought 
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that  this  should  give  a  better  intitial  presentation  of  the 
projections  than  any  other  choice  of  values.   However,  it 
should  be  pointed  out  that  the  decision  on  the  particular  set 
of  values  to  insert  is  left  up  to  the  user. 

Polaroid  prints  of  the  AGT  plot  of  the  projections  have 
been  included  with  the  computer  printout.   The  close  resemblance 
between  the  printout  and  the  Polaroid  pictures  is  quite 
apparent.   Of  course,  the  Polaroid  picture  is  a  true  repre- 
sentation of  the  projections.   Any  discrepancies  in  the  plot 
from  the  computer  printout  would  be  due  to  the  restrictions  in 
setting  up  the  matrix  array.   Adjustments  of  the  vertical  and 
horizontal  controls  were  necessary  to  center  the  plot  on  the 
screen  so  a  Polaroid  print  may  be  distorted  but  the  relative 
distance  between  points  has  been  maintained.   The  elaboration 
and  the  conclusion  of  each  case  will  now  be  presented. 

CASE  1 


The  first  case  involved  the  coordinates  of  a  three- 
dimensional  rectangular  pyramid.   The  points  represent  the 
outline  of  this  figure.   The  data  set  consisted  of  twenty-one 
vectors.   Each  consecutive  iteration  will  show  the  rotation 
of  the  geometrical  shape.   The  change  in  the  A  and  B  vectors 
merely  reoriented  the  perspective  of  the  figure.   Also  the 
shape  of  the  figure  remains  true  with  each  iteration.   This 
example  was  used  in  evaluating  the  successful  operation  of  the 
orthonormal  projections. 
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Figures  IV-1-1  through  IV-1-5  in  Appendix  V  show  a  five- 
step  sequence  in  the  application  of  the  orthonormal  projection 
concept.   Since  this  case  was  done  only  as  a  test  of  the 
algorithm,  no  printout  of  the  appropriate  terms  is  given.   The 
outline  of  this  geometrical  shape  was  drawn  in  by  hand  in  the 
figures.   The  Polaroid  picture  of  the  AGT  display,  Figure 
IV-1-5P  in  Appendix  V,  is  included  to  compare  with  the  print- 
out given  in  Figure  IV-1-5. 

CASE  2 

Two  three-dimensional  cubes  were  generated  for  analysis 
using  the  pattern  classification  technique.   The  X  and  Y 
subsets  comprised  twenty-tuples  depicting  the  outline  of  the 
two  geometrical  figures.   The  Y  subset  has  its  closest  vertex 
to  the  origin  at  the  point  (3,2,5).   Each  of  its  sides  has  a 
length  of  two  units.   The  other  cube  (X  subset)  had  its  closest 
vertex  at  (7,6,9)  and  also  two  units  on  a  side.   The  two  cubes 
were  made  disjoint.   This  made  it  feasible  to  specify  a  dis- 
criminant function  depicting  the  partitioned  subsets.   Again, 
this  was  another  test  in  evaluating  the  procedure  but  involved 
two  disjoint  figures. 

Figures  IV-2-1  through  IV-2-5  in  Appendix  V  show  a  six  step 
sequence  in  the  application  of  the  procedure  on  the  two  cubes. 
Separation  was  maintained  throughout  the  sequence.   The  values 
of  the  A  and  B  vectors  have  been  indicated  below  the  figures. 
The  initial  figure  used  the  values  from  the  table  in  Appendix 
IV. 
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A  possible  discriminant  function  in  Figure  I V -2-5  can 
be  specified  as 

f(Z)  =  A-(Z  -  DAV)/VMAX 

where  A  has  the  values  printed  from  Figure  IV-2-5.   If 
f(Z)  >  0,  Z  is  a  member  of  the  subset  of  Y's;  if  f(Z)  <  0, 
Z  is  an  entity  in  the  X's.   These  results  can  be  checked  by 
the  following  example.   The  printout  indicated  the  values 
for  VMAX  and  DAV ; 

DAV(l)  =  6.000        DAY(2)  =  5.000 
DAV(3)  =  8.000       VMAX    =  4.333 

Assume  that  Z  has  the  coordinates  of  the  point  Y(l)  in 
Appendix  III.   The  calculations  indicate  the  number  for  f(Z) 
is  +  0.662.   The  computation  has  shown  that  Z  fs ,  in  fact,  a 
member  of  the  subset  of  Y's.   Also  any  other  point  from  the 
Y's  will  satisfy  the  condition  that  f(Z)  will  be  positive. 
One  of  the  two  cubes  in  Figures  IV-2-4  and  IV-2-5  have 
been  drawn  in  by  hand.   Also  a  figure  of  the  Polaroid  print 
is  included  to  compare  with  Figure  IV-2-5. 

CASE  3 

For  the  third  application  using  the  new  separation  algo- 
rithm, the  dimensional  space  of  the  cube  was  increased  from 
three  to  six.   Seventy-tuples  were  generated  for  each  subset, 
X  and  Y.   The  unit  cubes  were  generated  on  the  IBM  360  computer 
as  uniform  distributed  random  numbers.   A  point  from  the  Y 
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subset  would  be  denoted  by  Y(I)  =  (Y(I,1),  Y(I,2),  Y(I,3), 
Y(I,4),  Y(I,5),  Y( I  ,6) )  where  I  =  1.....N  and  1  <  Y(I.K) 
<  2,  K  =  1,...,6;  and  for  the  X  subset  X(J)  =  (X(J,1),  X(J,2), 
X(J,3),  X(J,4),  X(J,5),  X(J,6))  where  J  =  1,...,M  and  0  < 
X(J,K)  £  1,  K  =  1,...,6.   Again,  the  two  cubes  were  generated 
so  the  two  samples  were  disjoint.   First,  if  one  of  the  itera- 
tions showed  a  combination  of  "X's"  and  "O's"  in  an  area  of 
the  plane,  it  would  be  easier  to  separate  the  two  groups  on 
subsequent  steps.   Secondly,  once  separati on  between  the  two 
groups  was  attained,  a  discriminant  function  could  be  defined. 

The  printout  of  the  graphical  display  is  shown  in  Figures 
IV-3-1  through  I V -3-6  in  Appendix  V.   The  degree  of  complexity 
for  the  two  samples  was  increased  from  that  of  Case  2,  but 
a  partition  of  the  two  groups  was  still  achieved.   Since  this 
case  is  so  similar  to  Case  2,  the  discriminant  function  will 
not  be  determined.   Consequently,  the  A  and  B  vectors  will 
also  be  omitted.   In  Figure  I V  -  3-2 ,  the  two  cubes  shown  by 
their  projections  can  not  be  separated  because  of  the  mixture 
of  "X's"  and  "O's"  around  the  center  of  the  graph.   It  was 
easy,  however,  to  separate  the  two  as  shown  in  the  next  step, 
Figure  IV-3-3.   The  similarity  between  the  printout  and  the 
Polaroid  snapshot  is  shown  in  Figures  IV-3-5  and  IV-3-5P. 

CASE  4 

The  previous  two  cases  were  applications  of  the  parti- 
tioning of  geometrical  shapes  which  were  separable  by  a  hyperplane 
A  sphere  and  a  shell  in  four-space  were  created  for  the  purpose 
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of  making  the  division  of  two  groups  more  complex.   Twenty- 
tuples  were  generated  for  the  X  subset  and  forty-tuples  for 
the  Y  subset.   A  typical  position  of  a  vector  in  the  Y  set  is 

Y(I)  =  (Y(I,1),  Y(1, 2),  Y(I,3),  Y(I,4)),  where  I  =  1,...,N. 

2         2 

The  components  must  satisfy  the  requirement  2  <  Y(I,1)   +  Y(I,2) 

+  Y(I,3)2  +  Y(I,4)2  <  4.   Similarly,  an  X  point  must  fulfill 
the  condition  that  X(J,1)2  +  X(J,2)2  +  X(J,3)2  +  X(J,4)2  <  1. 
The  two  subsets  are  again  disjoint  and  they  form  a  sphere 
surrounded  by  a  shell  separated  by  one  unit  in  radial  distance. 
Two  different  trials  on  the  sample  set  were  conducted.   The 
first  run  was  finally  terminated  after  twenty  iterations  without 
achieving  successful  separation.   It  should  be  noted  that  the 
arrangement  of  the  two  groups  made  it  impossible  for  an  appro- 
priate partition.   Of  course,  this  is  assuming  that  the  shell 
has  the  proper  density  of  points  and  no  voids  exist  in  it. 
However,  if  the  procedure  is  modified,  separation  can  be 
feasible.   The  second  run  was  carried  out  by  including  the 
square  of  each  component  in  each  entity.   The  dimension  of  the 
sample  set  becomes  twice  the  original.   To  clarify  this  last 

remark,  if  there  was  an  arbitrary  point,  Z  =  (Z,  ,Zp  ,Z-  ,Z . )  in 

2 
the  original  set,  it  will  now  become  Z  =  (Z,  ,Zp  ,Z3  ,Z .  ,Z,  , 

2    2    2 

Z  «  ,Z3  ,  Z .  ).   The  second  trial  was  conducted  using  this  eight 

component  vector  on  both  sets. 

The  first  trial  on  case  four  appears  in  Figures  IV-4-1 
through  IV-4-5  in  Appendix  V.   These  figures  confirm  that  the 
separation  was  never  attained.   Only  the  first  six  steps  have 
been  included  in  the  figures.   The  remaining  steps  were  omitted 
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since  they  did  not  improve  the  pattern  of  the  projections. 
Thus,  a  discriminant  function  can  not  be  specified.   Conse- 
quently, the  printout  of  the  A  and  B  vectors  will  be  eliminated 
from  the  discussion. 

The  situation  improves  in  the  second  run  using  the  eight 
component  tuples.   This  is  depicted  in  Figures  IV-4-6  through 
IV-4-10  in  Appendix  V.   Separation  was  attained  in  a  mere  three 
steps  of  the  algorithm.   This  modification  of  the  procedure 
can  be  attempted  on  samples  that  can  not  be  successfully 
partitioned  using  the  basic  method.   The  one  subset  may  be 
enclosed  inside  a  large  sphere-like  shape  including  the  other 
subset.   This  is  one  approach  to  obtaining  a  division  of  the 
two  groups.   The  components  of  the  A  and  B  vectors  in  the 
second  run  are  printed  on  the  bottom  of  each  of  the  figures. 

CASE  5 

This  sample  set  represents  the  first  of  three  "real  world" 
cases.   The  data  set  for  this  case  portrays  multiple  measure- 
ments on  three  species  of  Iris  flowers.   The  data  was  originally 
used  by  Ronald  A.  Fisher  in  his  classic  paper  on  discriminant 
analysis  [3],   It  has  become  widely  used  in  illustrating  cluster 
analysis  procedures.   The  four  measurements  collected  are 
(1  )  sepal  length,  (2)  sepal  width,  (3)  petal  length,  and  (4) 
petal  width.   The  sample  set  has  been  reduced  from  the  original 
150  to  75  samples  for  this  particular  analysis. 
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This  particular  sample  set  has  the  characteristics  for 
either  a  cluster  analysis  or  pattern  classification  problem. 
The  author  will  approach  the  investigation  using  the  pattern 
classification  technique.   It  is  quite  possible  from  looking 
at  the  values  of  the  three  species  that  the  one,  Iris  setosa, 
will  have  projections  in  an  area  of  the  plane  by  itself.   The 
projections  of  the  other  two  species  appear  to  be  nearly  the 
same.   Hence,  it  may  be  somewhat  difficult  to  partition  these 
latter  two  species.   Therefore,  the  observations  from  Iris 
versicolor  and  Iris  virginica  will  be  divided  into  two  different 
subsets.   The  Y  subset  will  be  composed  of  Iris  setosa  and 
Iris  versicolor;  while  the  X  subset  of  Iris  virginica.   Each 
species  of  flowers  was  composed  of  twenty-five  vectors  in 
four  dimensional  space. 

The  complexities  in  measurements  is  apparent  in  the  parti- 
tioning of  these  three  groups.   Outliers  may  exist  which 
indicate  an  observation  is  far  removed  from  the  main  body  of 
the  data.   They  result  from  the  fact  that  the  sample  under 
investigation  may  contain  a  certain  proportion  of  extreme 
cases  or  from  a  population  other  than  the  sample  from  which 
the  rest  of  the  data  came.   An  observation  from  the  Iris 
versicolor  species  presents  itself  as  a  possible  outlier  in 
the  analysis  of  the  graphic  plot.   This  point  is  Y(28)  in 
Appendix  III  which  coincides  with  the  projection  of  X(24). 
The  two  appear  as  one  point  in  the  last  two  figures  as  a  "J" 
symobl .   The  point  Y(28)  could  not  be  separated  from  the  group 
of  X's  as  it  is  projected  to  the  left  of  the  line  denoting 
the  discriminant  function. 
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Figures  IV-5-1  through  IV-5-6  in  Appendix  V  show  the 
twelve  step  sequence  in  the  analysis  of  the  three  Iris 
flowers.   Even  though  a  discriminant  function  is  generally 
determined  in  partitioning  two  groups,  two  linear  functions 
will  be  specified  which  divide  the  three  types  of  flowers. 
This  will  be  done  using  Figure  IV-5-6.  The  first  function 
will  separate  the-  X  subset  from  the  Y  subset.   Then  the 
second  will  divide  the  two  Iris  species  comprising  the  Y 
subset.   First,  the  values  of  the  A  and  B  vectors  from  the 
bottom  of  the  figure  must  be  recorded. 

Now,  a  linear  discriminant  function  has  the  form  given  in 
Secti  on  II  as ; 

-m    d(A(Z  -  DAV))     e(B(Z  -  DAV))  .  f 
Tu;  ~      VMAX  VMAX 

where  d,  e,  and  f  are  scalars.   The  f  scalar  is  the  negative 
sign  of  the  ordinate  intercept  while  the  term  d  has  the 
negative  sign  of  the  slope  of. the  line.   Finally,  the  value 
of  the  coefficient  e  will  be  one.   Using  the  line  in  Figure 
IV-5-6  separating  the  X  and  Y  groups,  an  appropriate  dis- 
criminant function  has  the  form: 

f(7)    -    6(A(Z  -  DAV))    B(Z  -  DAV)    2  25 
T{L)  VMAX  VMAX        ' 

The  magnitude  for  DAV  and  VMAX  from  this  figure  are: 

DAV(l)  =  5.872  DAV(2)  =  3.061 

DAV(3)  =  3.804  DAV(4)  =  1.212 
VMAX    =  3.033 
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Using   these    figures,    the    value    for   f(Z)    is    -    0.054   when   the 
components    of   X(24)    from   Appendix    III    are    substituted   into    Z. 
Also   carrying   out    the   same    computations,    the    result    for    f(Z) 
is    +   0.252    for   the    point   Y(26).      These    calculations    have 
merely   emphasized    that   a    positive    value    for    f(Z)    implies    Z 
is    a   member   of   the   Y    subset.      If   f(Z)    is    negative,    then   Z    is 
an   entity   of   the    X   subset. 

Now,    a    description    of   the   line    splitting   the   Y    subset 
must   be    resolved.      The    form   of   the    discriminant    function   has 
been    determined   as: 

f(U)    =   A'(U    -    DAV)    _    o    370 
UUj  VMAX  U.J/U 

where    the    number   0.370    is    the    value    of   the    intercept   on    the 
abscissa    axis.      If   the    coordinates    of   Y(38)    are    substituted 
for   the   arbitrary   point    U,    the    result   of   f(U)    is    -    0.030. 
This    calculation    indicates    that    any    point    from   the    set   of 
Y's,    where    f(U)    <   0,    represents    an    entity    from   the    Iris 
versicolor   species.      Any    U   where    f(U)    >   0    corresponds    to   an 
observation    from   the    Iris    setosa. 

CASE    6 

The  sample  set  for  this  analysis  problem  is  quantitative 
measurements  of  morphologic  parameters  of  specific  fossils. 
The  data  set  was  extracted  from  Wright  and  Switzer's  paper 
titled,  "Numerical  Classification  Applied  to  Certain  Jamaican 
Eocene  Nummul i tids  "  [4].   Eighty-two  eight  component  vectors 
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were  used  in  the  investigation.   Six  of  the  original  observa- 
tions had  missing  parameters  and  were  therefore  omitted  from 
the  sample  set. 

Again  like  Case  5,  this  sample  will  be  divided  into  two 
groups  for  a  pattern  classification  problem.   The  author  feels 
the  separation  will  be  easier  to  attain  using  two  subsets. 
It  should  be  noted  that  Lauzon  investigated  these  specimens 
using  cluster  analysis  procedures  [1].   The  Y  subset  was 
initially  composed  of  the  first  thirty-eight  observations. 
It  is  comprised  of  specimens  from  the  Opercul inoi des  wilcoxi 
species.   This  species  will  be  denoted  by  Group  A.   The  X 
subset  was  composed  of  the  other  three  species:   Opercul i noi des 
jennyi  ,  Opercul i noi des  cushmani  ,  and  Nummulites  striatoreti- 
culatus.   These  three  types  of  Nummulitids  will  be  called 
Groups  B,  C,  and  D,  respectively,   Using  these  subsets,  it 
was  extremely  difficult  to  pick  out  the  projections  from 
Group  B  and  those  from  Group  C.   In  order  to  eliminate  this 
confusion,  the  two  subsets  were  rearranged.   Observations 
from  Groups  A  and  B  were  combined  to  form  the  Y  subset; 
Groups  C  and  D  joined  to  make  up  the  X  subset.   The  second 
run  was  conducted  using  the  above  groupings.   Also  to  make 
the  investigation  similar  to  the  discussion  in  Wright  and 
Switzer's  paper,  only  the  first  six  parameters  were  used. 
The  results  from  the  analysis  were  quite  remarkable.   Even 
on  the  first  step,  three  different  groups  were  apparent; 
Group  A,  Groups  B  and  C,  and  Group  D.   In  Wright  and  Switzer's 
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investigation,  ninety-six  iterations  were  conducted  to  deter- 
mine the  coefficients  of  a  linear  discriminant  function.   They 
were  able  to  find  two  different  groupings,  Groups  A  and  D  and 
Groups  B  and  C.   The  same  procedure  using  96  iterations  was 
repeated  to  find  another  linear  function  dividing  Group  A  and 
D.   At  last,  a  linear  function  splitting  Group  B  and  C  was 
determined.   Then-  the  two  researchers  utilized  these  two  final 
divisions  as  input  to  a  standard  discriminant  analysis  program, 
The  output  resulted  in  two  canonical  variables.   The  plot  of 
these  two  variables  revealed  three  groupings:  Group  A,  Groups 
B  and  C,  and  Group  D.   The  two  concluded  that  the  specimens 
from  Groups  B  and  C  could  not  be  partitioned.   The  results 
in  this  thesis  are  in  agreement  with  these  findings.   The 
procedure  used  by  this  author  appears  to  be  superior  in  two 
ways.   First,  there  are  substantially  fewer  computations, 
hence  the  computer  time  to  accomplish  the  analysis  is  consider 
ably  less.   Secondly,  the  plot  of  the  results  is  readily 
available  from  the  printout  of  the  computer.   It  should  also 
be  noted  that  G.  Lauzon's  investigation  of  this  sample  set 
showed  only  two  groupings. 

Figures  IV-6-1  through  IV-6-3  in  Appendix  V  show  the 
graphical  display  obtained  on  the  first  run.   These  figures 
depict  the  eight-dimensional  analysis  of  the  fossil  measure- 
ments.  In  Figure  IV-6-2,  the  three  groups  have  been  labeled. 
There  is  no  visual  method  to  distinguish  between  the  B  and  C 
groups.   The  displays  of  the  second  run  are  given  in  Figures 
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IV-6-4   through    IV-6-7.      The    values    of   the   A   and   B  components 
of   these    figures    are   given   below   each. 

The    two    discriminant    functions    separating    the  three 

groupings   will    be   specified   using    Figure    IV-6-7.  First,    the 

tabulation   of   the   numbers    for   VMAX    and   DAV   are    as  follows: 

DAV(l)  =  163.220  DAV(2)  =  41 .43 

DAV(3)  =  9.390  DAV(4)  =  25.012 

DAV(5)  =  88.866  DAV(6)  =  677.817 

VMAX    .  =  587.954 

A   description   of   the    line   splitting    Groups    A   and   D   must   be 
determined.      The    form   of   the    discriminant   function   has    been 
resolved   as 

f(7\    -    B*(Z    -    DAV)        n    ,-,. 

If  f(Z)  >  0,  then  the  arbitrary  point  Z  is  a  member  of  Group 

A,  B,  or  C.   On  the  other  hand,  if  f(Z)  <  0,  then  Z  is  an 

entity  of  Group  D.   The  equation  of  the  discriminant  function 
describing  the  upper  line  is  given  as: 

-nn      0.6(A-(U  -  DAV))    B-(U  -  DAV) 
nu;    "         VMAX  VMAX 

U  is  an  observation  from  Group  A  or  D  if  f(U)  >  0.   Likewise, 
U  is  a  vector  from  Group  B  or  C  if  the  result  of  f(U)  <  0. 
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CASE  7 

The  final  analysis  was  conducted  on  a  thi rteen-dimensi onal 
pattern  classification  problem.   Coronary  patients  from 
Letterman  Hospital  comprised  the  sixty-tuples  making  up  the 
Y  subset  [1].   The  other  subset  was  composed  of  fifty  Army 
personnel  at  Silas  B.  Hays  Hospital  who  had  undergone  annual 
physical  examinations.   Those  patients  from  Letterman  had 
been  diagnosed  as  having  some  type  of  coronary  heart  disease. 
The  X  subset  consisted  of  thirteen  features  from  a  random 
sample  of  the  Army  personnel.   There  were  no  indications  from 
the  latter  group  to  either  classify  or  not  as  having  heart 
disease.   If  successful  separation  is  obtained  and  a  discrimi- 
nant function  is  computed,  then  it  could  be  used  to  classify 
new  patients. 

A  few  of  the  components  in  each  thi rteen-tupl e  represented 
measurable  entities  such  as  weight  and  age.   Other  character- 
istics such  as  race  and  smoking  behavior  were  assigned 
arbitrary  values  rather  than  numbers  to  signify  a  measurement. 
The  analysis  will  show  that  these  classifiers  are  not  really 
appropriate  to  use.   The  actual  measurements  exhibit  at  least 
the  characteristics  of  an  interval  scale.   This  is  a  minimum 
requirement  in  computing  a  meaningful  average  value  of  a 
parameter.   Those  components  which  had  an  arbitrary  value 
assigned  do  not  meet  this  criteria  for  an  interval  scale.   In 
fact,  they  have  the  characteristics  of  an  ordinal  scale. 
Therefore,  when  the  average  values  of  the  components,  DAV, 
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are  computed  in  the  program,  it  is  inappropriate  to  use  these 
classifiers  in  the  analysis. 

The  first  run  was  conducted  on  the  thi rteen-dimensi onal 
sample  set.   The  displays  of  the  projections  are  given  in 
Figures  IV-7-1  through  IV-7-6.   The  last  figure  had  the  best 
presentation  with  respect  to  partitioning  the  two  groups.   A 
line  dividing  the"X's"  and  "O's"  was  drawn  in  the  figure. 
If  a  discriminant  function  had  been  specified  for  this  line, 
thirty-six  of  the  fifty  X  points  and  forty-nine  of  the  sixty 
Y  points  would  have  been  properly  classified.   This  result 
means  that  approximately  25%  of  the  total  sample  have  been 
projected  out  of  their  respective  regions.   This  last  state- 
ment does  not  seem  satisfactory  for  any  prediction  model  or 
separation . 

The  values  of  the  A  and  B  vectors  used  by  G.  Lauzon  of 
this  data  indicate  a  definite  pattern.   For  those  assigned 
values,  the  number  in  both  vectors  for  each  component  are  less 
than  or  equal  to  the  absolute  value  of  two.   On  the  other 
hand,  one  of  these  two  numbers  is  wery   much  greater  than  2 
for  every  measured  component.   For  this  reason,  the  components 
of  assigned  values  will  be  deleted  in  the  second  run.   In 
addition,  this  algorithm  is  more  appropriate  for  separation 
of  data  where  the  sample  set  is  on  an  interval  scale. 

In  the  second  run,  an  arbitrary  observation  Z  will  be 
denoted  by 

Z(H)  =  (Z(H,2),  Z(H,3),  Z(H,4),  Z(H,12),  Z(H,13)) 
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where   H    =    1,...,M   for   an    X   and  H    =    1,...,N   for   a   Y.      Figures 
IV-7-7    through    IV-7-12    show    the    printout    of   the    five-dimensional 
plot.      The    figure   with    the    best    partitioning    is    IV-7-12.      This 
time    the    discriminant    function   would    result    in    a    smaller 
probability   of  error.      There    are    fifty-four   Y    projections    in 
the    left    region    and   thirty-six    X's    to    the    right.      Now,    the 
result   is    less    than   80%    of   the    projections    fell    into    their 
respective    areas.      The    author   feels    that    this    second    run   has 
also    led    to    an    unsatisfactory   separation   of   the    two    groups. 
The    reason    for   this    conclusion   is    the    high    density    of   points 
near   the    dividing    line.      In    addition,    the    distance    is    rather 
insignificant   for   the    different   points    of   the    two   subsets 
across    this    dividing    line.      Therefore,    the    author   feels    that 
an   attempt    to   specify   a    discriminant   function    for   this    data 
would    be    futile. 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

It  can  be  concluded  that  the  improvements  on  the  algorithm 
have  resulted  in  achieving  greater  success  than  the  previous 
application.   First  of  all,  the  revised  algorithm  provided  a 
partitioning  of  the  generated  sample  sets  in  the  first  four 
cases.   The  new  separation  algorithm  was  also  applied  to  three 
"real  world"  applications.   The  first  investigation  obtained 
a  satisfactory  partitioning  of  the  three  Iris  species.   The 
three  groups  were  divided  by  specifying  two  discriminant 
functions.   The  second  sample  set  involved  a  more  complex 
analysis  problem.   The  observations  consisted  of  an  eight- 
dimensional  array.   In  addition,  it  had  been  hypothesized 
that  there  existed  four  species  in  the  data.   The  display  for 
this  analysis  indicated  the  presence  of  only  three  species. 
The  same  conclusion  was  also  reached  by  Wright  and  Switzer  [4] 
in  their  approach  to  this  discriminant  analysis  problem. 
However,  their  computations  in  reaching  this  result  were  much 
more  involved  than  the  algorithm  in  this  paper.   The  investi- 
gation by  G.  Lauzon  showed  the  presence  of  only  two  such 
separate  groupings.   From  this,  the  improved  algorithm  using 
the  orthonormal  vectors  appears  to  achieve  a  more  definite 
partitioning.   Finally,  the  employment  of  this  method  on  the 
thi rteen-dimensi onal  medical  data  resulted  in  unsatisfactory 
separation.   The  rejection  of  a  partition  may  be  due  to  the 
use  of  inappropriate  parameters  in  analyzing  coronary  disease. 
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All  separations  are  really  approximations.   The  final 
classifications  must  be  subjected  to  a  variety  of  sample  sets 
for  a  confirmation  of  the  discriminant  function.   Even  after 
this  has  been  accomplished,  the  presence  of  outliers  will 
result  in  some  finite  value  for  the  probability  of  error.   The 
description  by  the  discriminant  function  is  only  the  optimal 
division  of  the  groupings. 

The  use  of  the  values  from  the  table  in  Appendix  IV  should 
be  mentioned.   All  investigations  in  this  paper  were  conducted 
choosing  the  initial  orthonormal  vectors  as  listed  in  the  table 
(Appendix  IV).   However,  if  different  values  had  been  sub- 
stituted in  the  components,  the  final  separation  may  be  better. 
In  the  future,  the  operator  should  utilize  a  variety  of  unit 
orthonormal  vectors  in  the  attempt  of  achieving  a  successful 
parti  ti  on . 

This  algorithm  has  provided  a  relatively  easy  and  simple 
approach  to  the  classification  problem.   It  is  possible  that 
the  merit  of  this  method  could  be  incorporated  into  a  portable 
unit.   A  computer  module  designed  for  the  appropriate  computa- 
tions would  provide  a  means  for  obtaining  the  required 
projections.   A  television  cathode  ray  tube  could  act  as  the 
display  system.   An  interface  system  could  be  developed  to 
handle  the  input/output  operations.   The  outstanding  feature 
using  a  color  t.v.  set  is  the  ability  to  display  a  variety 
of  colors.   The  various  subsets  of  the  sample  set  could  be 
linked  to  a  respective  hue.   The  combination  of  these  two 
systems  would  make  up  an  ideal  graphics/computer  system  in 
which  the  algorithm  would  be  a  major  feature. 
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The  development  of  this  algorithm  could  be  applied  in  an 
unlimited  variety  of  fields.   The  analysis  of  medical  and 
taxonomy  samples  has  already  been  mentioned.   The  method  is 
pertinent  to  any  field  in  which  there  is  mutivariate  data  to 
be  separated.   For  instance,  measurements  of  meteorological 
parameters  may  be  subjected  to  this  method  of  analysis.   The 
motivation  of  the  investigation  could  be  to  determine  the 
existence  of  conditions  pertaining  to  the  development  of  a 
tropical  storm  and  those  present  during  the  infancy  of  a 
hurricane.   This  algorithm  may  be  able  to-  distinguish  the  two 
different  groupings.   If  this  can  be  accomplished,  then  the 
model  can  be  used  as  a  basis  in  predicting  future  storm 
acti vi  ty . 

The  final  comments  are  directed  to  future  changes  in  the 

computer  routine.   First,  the  A  and  B  vectors  should  be 

punched  on  computer  cards  and  read  in  with  the  data  deck. 

This  would  eliminate  the  possibility  of  typing  in  an  incorrect 

component.   Secondly,  a  video  presentation  of  the  number  of 

the  selected  point  would  be  most  helpful  to  the  operator. 

In  addition,  the  printout  of  the  plot  should  indicate  the 

selected  point  by  a  different  set  of  letters  than  the  present 

"X-0"  scheme.   These  two  modifications  would  make  it  possible 

to  repeat  the  same  iterations  on  subsequent  runs.   Lastly, 

it  would  be  beneficial  to  the  operator  if  a  series  of  sample 

sets  could  be  read.   This  could  be  arranged  by  adding  the  use 

of  another  function  switch.   The  operator  would  now  be  able  to 

conduct  a  series  of  different  investigations  without  leaving 

the  graphics  console. 
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APPENDIX  II.   PROGRAM  OPERATING  INSTRUCTIONS 

There  exists  two  methods  by  which  this  program  can  be 
executed  on  the  XDS  9300.   The  first  and  fastest  method  is 
utilizing  a  magnetic  tape.   The  memory  dump  of  the  compila- 
tion from  the  digital  computer  has  been  recorded  onto  a 
magnetic  tape.   The  elapsed  time  until  execution  commences 
would  be  approximately  five  minutes  using  this  tape.   The 
other  method  is  simply  reading  the  Fortran  program  from 
the  computer  cards.   It  would  take  roughly  thirty  minutes 
until  the  program  is  ready  to  execute  using  the  cards.   This 
time  factor  is  very  crucial  if  more  than  one  set  of  data 
cards  must  be  read. 

In  executing  the  program  from  the  magnetic  tape,  the 
operator  may  consult  the  publication  titled  "Electrical 
Engineering  Computer  Laboratory"  by  R.  D.  Delaura  [2].   A 
brief  guide  in  using  the  magnetic  tape  mechanism  applicable 
to  the  object  language  filed  on  this  tape  follows: 

(1)  Load  the  magnetic  tape 

(2)  Unit  select  to  position  "1" 

(3)  Automatic/manual  switch  to  "MANUAL" 

(4)  Depress  forward  drive  (tape  will  stop  at 
the  LOAD  POINT) 

(5)  Automatic/manual  switch  to  "AUTOMATIC" 

The  magnetic  tape  is  now  ready  for  instructions  from  the 
XDS  9300.   The  operator  must  incorporate  five  control  cards 
into  the  data  deck.   These  control  cards  are  instructions 
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to  the  various  systems  used  in  the  execution  of  this  program. 
The  deck  of  cards  will  be  as  follows: 

(1)  B00T 

(2)  J0B 

(3)  AGT 

(4)  ASSIGN  XI  =  MT1A 

(5)  RERUN 

(6)  (DATA  DECK) 

The  deck  is  now  ready  to  load  into  the  card  reader  and 
execution  can  commence. 

The  following  instructions  will  be  applicable  to  the 
use  of  computer  cards.   The  list  of  control  cards,  that 
must  precede  the  source  program  printed  in  Appendix  I, 
are  given  below: 

(1)  B00T 

(2)  $  DATA 

(3)  $PATCH 

(4)  007067  001  06711 

(5)  $END 

(6)  J0B 

(7)  AGT 

(8)  F0RTRAN  LS,  G0 

The  data  deck  is  placed  at  the  end  of  the  Fortran  program. 
The  deck  of  cards  is  now  ready  to  load  and  execution  can 
commence . 

In  preparing  the  computer  program  using  either  method, 
four  cards  must  be  punched  and  inserted  into  the  data  cards. 
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The  first  three  cards  can  be  used  for  the  purpose  of  verbally 
describing  the  data  to  be  investigated.   This  information 
can  be  punched  using  any  part  or  all  of  the  eighty  column 
card.   The  operator  must  insure  in  insert  these  three  cards  even 
though  one  or  more  cards  may  be  blank.   The  fourth  card  is  the 
instruction  by  which  the  computer  reads  the  data  from  the 
punched  cards.   For  illustrative  purposes,  The  statement  might 
read,  "  (4F7,2 ,3X,3F9 . 2)  . "   This  indicates  that  a  total  of 
seven  parameters  will  be  read  from  each  card.   The  first 
parenthesis  must  be  punched  in  column  seven  of  the  card. 
The  actual  data  cards  are  inserted  behind  these  four  cards. 
The  four  cards  plus  the  cards  of  data  will  be  considered 
the  data  deck. 

A  guide  for  operating  the  XDS  9300  computer  and  the  AGT 
console  using  either  a  magnetic  tape  or  cards  will  be  given. 
First  though,  assume  that  there  are  M  number  of  X  points 
and  N  number  of  Y  points.   The  value  of  M  and  N  must  be 
less  than  or  equal  to  eighty.   In  addition,  all  observations 
in  both  sets  constitute  an  L-dimensi onal  vector.   The  number 
of  L  must  be  less  than  or  equal  to  thirteen.   Now,  the 
instructions  will  be  presented. 

(1)  Use  the  "POWER  ON"  procedure  for  the  XDS  9300 
and  the  AGT  graphic  systems. 

(2)  Executing  the  program 

(a)  At  the  XDS  9300  control  panel,  depress  "SENSE 
2"  if  using  cards.   Do  not  depress  if  using  the 
tape.   Then  depress  "IDLE",  "CLEAR"  and  "CLEAR 
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FLAGS"  simultaneously,  "RESET",  "RUN",  and 
"CARDS."  The  card  reader  will  commence  reading 
the  punched  cards.   The  line  printer  and 
teletypewriter  (TTY)  will  simultaneously  print 
"JOB."  The  computer  will  now  start  compiling 
the  program. 

(b)  At  the  time  execution  of  the  program  starts,  the 
TTY  will  print  the  message:   "TYPE  IDEV  =  2*  AND 
A  CARRIAGE  RETURN  IF  USING  AGT2."   The  input 
light  will  illuminate  to  alert  the  operator. 
The  operator  must  type  the  number  of  X's  and  Y's 
plus  the  dimension  of  the  two  subsets  of  data. 
In  addition,  the  operator  has  the  option  of 
printing  all  WRITE  statements  prior  to  statement 
#57  in  the  computer  program.   This  is  done  by 
typing  "NPFLG  =  1"  if  a  printout  is  desired; 
NPFLG  =  0"  for  no  printout.   As  an  example, 
suppose  there  are  15  X's  and  20  Y's  in  10 
dimensions.   If  the  operator  is  using  AGT  2  and 
requests  no  printout,  type  "M  =  15,  N  =  20,  L  =  10, 
NPFLG  =  0,  IDEV  =  2,  *"  and  a  carriage  return  at 
the  computer  TTY.   If  using  AGT  1,  simply  delete 
"IDEV  =  2."   The  card  reader  will  begin  the 
process  of  reading  the  card. 

(c)  The  word  "NAMELIST"  and  a  blinking  cursor  will 
be  displayed  on  the  bottom  left  corner  of  the 

AGT  screen.   This  will  occur  when  computer  reaches 
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statement  #57  in  the  program.   The  operator 
must  type  in  the  components  of  the  A  and  B 
vectors.   The  values  from  Appendix  IV  may  be 
used  or  the  operator  may  use  any  two  orthonormal 
unit  vectors.   Thirteen  numbers  for  each  vector 
must  be  inputted.   If  the  sample  is  10-di mensi onal , 
zero.es  will  be  inserted  for  the  last  three 
components.   A  representation  of  the  above  example 
would  look  like:   "A  =  .31628,  .31628,...,  .31628, 
0,0,0"  and  a  carriage  return,  "B  =  .31628,..., 
-.31628,0,0,0"  and  a  carriage  return,  "*"  and  a 
carriage  return.   The  XDS  9300  will  compute  the 
projections  and  relay  the  data  to  the  AGT  for 
display  on  the  screen.   The  "X's"  and  "O's"  denote 
the  original  "X's  and  Y's,  respectively. 
(3)  Interations  and  output 

(a)  To  change  the  orientation  of  the  plane,  the 
operator  chooses  either  the  "PICK  X"  or  "PICK 
Y"  function  switches.   The  location  of  these 
buttons  on  the  overlay  are  shown  in  Figure 
A-2-1.   A  video  message,  "TEXT  EDIT  REQUEST," 
will  appear  at  the  bottom  of  the  screen.   The 
operator  depresses  the  button  on  the  lightpen 
and  places  the  light  over  the  desired  symbol 
(point)  to  be  moved.   Upon  releasing  the  button, 
a  portion  of  the  letter  will  disappear  to 
signify  the  selected  point.   The  user  pushes 
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the  two  switches  "COMPLETE  CHOICE"  and  "END 
EDIT".  The  video  message  will  disappear.   The 
iteration  will  commence  after  one  of  the  move- 
ment buttons  (function  switches)  is  depressed: 
"RIGHT",  "LEFT",  "UP",  "DOWN."   The  9300 
recomputes  the  projections  from  the  inner  product 
with  the  orthonormal   vectors  and  displays  the  new 
results  on  the  screen. 
(b)  A  printout  of  the  new  projections  can  be  obtained 
b-y  simultaneously  depressing  the  "PRINT"  and  one 
of  the  movement  function  switches.   The  movement 
switch  must  be  released  first.   A  copy  of  the  new 
projections  on  the  plane  will  be  printed  at  the 
line  printer.   If  the  plot  of  the  computer 
projections  is  desired,  then  depress  the  "PLOT" 
f uncti  on  swi  tch . 
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Figure  A-2-1  .   FUNCTION  SWITCH  OVERLAY 
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APPENDIX  III.   DATA  LISTING 

This  appendix  presents  a  listing  of  the  data  sets  used 
in  the  seven  cases  discussed  in  Section  IV.   The  data  has 
been  arranged  so  that  a  row  of  numbers  represents  a  single 
vector.  For  example,  the  data  set  for  Case  one  is  presented 
in  three  columns  denoting  a  three-dimensional  entity.   The 
thirteen  columns  in  Case  Seven  represents  a  thirteen- 
dimensional  problem.   As  an  aid  for  reference  purposes,  the 
following  index  is  provided: 


NUMBER 

NUMBER 

PAGE 

OF  Y'S 

OF  X'S 

CASE  1 

21 

'  0 

CASE  2 

20 

20 

CASE  3 

70 

70 

CASE  4 

40 

20 

CASE  5 

50 

25 

CASE  6 

38 

44 

CASE  7 

60 

50 
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CASE      1       DATA 


Y(I,1> 


Yd, 2) 


Yd, 3) 


1 

1« 

0000 

1.0000 

1.0000 

2 

1< 

0000 

1.0000 

2.0000 

3 

2« 

5000 

1.0000 

1.0000 

h 

2. 

.5000 

1.0000 

2.0000 

5 

<H 

0000 

1.0000 

1.0000 

6 

4t 

.0000 

1.0000 

2.0000 

7 

1< 

.0000 

3.0000 

1.0000 

8 

1. 

.0000 

3.0000 

2.0000 

9 

2. 

.5000 

3.0000 

1.0000 

10 

2. 

.5000 

3.0000 

2.0000 

11 

4< 

.0000 

3.0000 

1.0000 

12 

4. 

.0000 

3.0000 

2.0000 

13 

1« 

.0000 

2.0000 

1.0000 

H 

1< 

•  0000 

2.0000 

2.0000 

15 

4. 

.0000 

2.0000 

1.0000 

16 

4, 

.0000 

2.0000 

2.0000 

17 

2« 

.5000 

2.0000 

3.3229 

18 

1 

►  7500 

1.5000 

2.6614 

19 

1 

.7500 

2.5000 

2.6614 

20 

3 

►  2500 

1.5000 

2.6614 

21 

3 

.2500 

2.5000 

2.6614 
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CASE      2      DATA 


Y(I#1) 


Y(N2) 


Y(I/3) 


1 

3.0000 

2.0000 

5.0000 

2 

4.0000 

2.0000 

5#0000 

3 

5.0000 

2.0000 

5,0000 

4 

3.0000 

4.0000 

5.0000 

5 

4.0000 

4.0000 

5,0000 

6 

5.0000 

4.0000 

5.0000 

7 

5.0000 

3.0000 

5.0000 

8 

3.0000 

3.0000 

5.0000 

9 

3.0000 

2.0000 

7.0000 

10 

4.0000 

2.0000 

7.0000 

11 

5.0000 

2.0000 

7.0000 

12 

3.0000 

3.0000 

7.0000 

13 

3.0000 

4.0000 

7.0000 

14 

4.0000 

4.0000 

7.0000 

15 

5.0000 

4.0000 

7.0000 

16 

5.0000 

3.0000 

7.0000 

17 

5.0000 

2.0000 

6.0000 

18 

5.0000 

4.0000 

6.0000 

19 

3.0000 

4.0000 

6.0000 

20 

3.0000 

2.0000 

6.0000 
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CASE      2       (C0NT) 


.  * 


XU/1) 


X(Ji2) 


X(J,3) 


1 

7.0000 

6.0000 

9.000 

a  . 

8.0000 

6.0000 

9.000 

3 

9.0000 

6.0000 

9.000 

f 

7.0000 

8.0000 

9.000 

5 

8.0000 

8.0000 

9.000 

6 

9.0000 

8.0000 

9.000 

7 

9.0000 

7.0000 

9.000 

8 

7.0000 

7.0000 

9.000 

9 

7.0000 

6.0000 

11.000 

10 

8.0000 

6.0000 

11.000 

11 

9.0000 

6.0000 

11.000 

12 

7.0000 

7.0000 

11.000 

13 

7-0000 

8.0000 

11.000 

14 

8.0000 

8.0000 

11.000 

15 

9.0000 

8.0000 

11.000 

16 

9.0000 

7.0000 

11.000 

17 

9.0000 

6.0000 

10.000 

18 

9.0000 

8.0000 

10.000 

19 

7.0000 

8*0000 

10.000 

20 

7.0000 

6.0000 

10.000 
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CASE      3      DATA 


YU/1) 


YC/2) 


Y(I/3) 


Y(I*4) 


Yd, 5) 


Y(Ii6) 


1.180              1 

1.718 

1.100             1 

.892 

1.617 

1.703 

1.597              1 

L.621 

1.526              1 

.358 

1.552 

1.746 

ltl97            : 

1.597 

1.540              1 

.590 

1.733 

1.484 

1.506              J 

1.561 

1.035              1 

.176 

1.892 

1.899 

1«588     .       : 

1.328 

1.639              1 

.625 

1.754 

1.684 

1.018              1 

1.094 

1.289              1 

.646 

1.284 

1.520 

1*517              J 

,.627 

1.317              1 

.812 

1.197 

1.137 

1.701              3 

..452 

1.342              1 

.452 

1.870 

1.917 

1.124              1 

,.239 

1.711              1 

.215 

1.546 

1.209 

1.883              J 

..120 

2.000              1 

.072 

1.344 

1.397 

1.649              1 

L.420 

1.998              1 

.861 

1.498 

1.810 

1.032            j 

L.365 

1.913              1 

.283 

1.460 

1.825 

1.929             : 

1.224 

1.880              1 

.747 

1.405 

1.080 

1.064              1 

L  .996 

1.541              1 

.837 

1.046 

1.039 

1.069              J 

1.291 

1.979             1 

.638 

1.170 

1.195 

1.092              2 

L.226 

1.174              1 

.806 

1.244 

1.802 

1.524              1 

L  .337 

1.919              1 

•  166 

1.532 

1.948 

1.176              \ 

L.690 

1.487              1 

.687 

1.422 

1.413 

1.585             : 

L  .648 

1.962              1 

.225 

1.378 

1.627 

1*632             : 

L.175 

1.615              1 

.571 

1.724 

1.216 

i.i7i          ; 

L.659 

1.643              1 

.252 

1.842 

1.958 

1.684          ; 

L  .932 

1.437              1 

.162 

1.378 

1.827 
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CASE      3       (C8NT) 
YCI/1)  Y(I/2)  Y(I/3) 


Y<I>4> 


Yd, 5) 


Y(I>6) 


23 

1«208              1 

.206 

1.544 

1.267 

1.497 

1.653 

2* 

1.968              1 

.342 

1.53J 

1.487 

1.153 

1.074 

25 

1.490              1 

..133 

1.575 

1.929 

1.084 

1.715 

26 

1.692              1 

..328 

1.899 

1.830 

1.106 

1.222 

27 

1.530              1 

,.909 

1.522 

1.227 

1.805 

1.239 

28 

1.003              I 

..158 

1.080 

1.096 

1.783 

1.462 

29 

1.344              1 

,.059 

1.792 

1.204 

1.515 

1.290 

30 

1.649              J 

..484 

1.656 

1.244 

1.874 

1.487 

31 

1.461              J 

,.181 

1.636 

1.397 

1.924 

1.152 

32 

1.701              J 

..317 

1.784 

1.409 

1.747 

1.171 

33 

1.745              3 

[.655 

1.551 

1.389 

1.368 

1.882 

34 

1.386              1 

1.039 

1.643 

1.236 

1.158 

1.472 

35 

1.757             i 

L  .371 

1.398 

1.550 

1.861 

1.213 

36 

1.122              J 

L  .396 

1.045 

1.647 

1.294 

1.450 

37 

1.289             : 

[.639 

1.797 

1.641 

1.825 

1.483 

38 

1.848             : 

1.173 

1.978 

1*774 

1.545 

1.693 

39 

1.519 

1.047 

1.679 

1.418 

1.915 

1.388 

40 

1.312 

L.741 

1.703 

1.528 

1.054 

1.775 

41 

1.777 

L  .019 

1.415 

1.232 

1.442 

1.839 

42 

1.504 

1.072 

1.080 

1.236 

1.507 

1.236 

43 

1.489 

1.230 

1.244 

1.393 

1.287 

1.450 

44 

1.599 

1.176 

1.485 

1.616 

1.422 

1.633 
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CASE      3      (C6NT) 
Yd,l)  Y(I/2)  Yd, 3) 


Yd, 4) 


Yd, 5) 


Yd, 6) 


1.997            : 

1.119 

1.400 

1.364 

1.144 

1.194 

1.239            : 

L.355 

1.864 

1.233 

1.052 

1.657 

1#842             : 

1.195 

1.941 

1.176 

1.052 

1.950 

1.724              ; 

1.840 

1.602 

1.833 

1.907 

1.756 

1.007              J 

1.528 

1.540 

1.648 

1.532 

1.372 

1.952     .        j 

L.902 

1.775 

1.218 

1.930 

1.344 

1.668            : 

1.392 

1.636 

1.575 

1.292 

1.937 

1.395              3 

1.263 

1.838 

1.016 

1.905 

1.538 

1.492              ) 

L.989 

1.514 

1.500 

1.607 

1.803 

1.629              J 

L  .333 

1.027 

1.040 

1.741 

1.282 

1.977              3 

L.646 

1.145 

1.948 

1.372 

1.585 

1.474              J 

[.988 

1.088 

1.578 

1.882 

1.557 

1.904              J 

..334 

1.104 

1.974 

1.088 

1.021 

1.184              J 

..770 

1.259 

1.116 

1.551 

1.654 

1.029             : 

..262 

1.112 

1.501 

1.344 

1.350 

1.582              1 

..701 

1.598 

1.656 

1.025 

1.647 

1.592              1 

,.496 

1.963 

1.020 

1.805 

1.703 

1.296             3 

.101 

1.369 

1.092 

1.789 

1.297 

1.059             1 

..146 

1.531 

1.058 

1.992 

1.649 

1.165             ] 

..174 

1.485 

1.724 

1.140 

1.005 

1.430             3 

,.907 

1.0S9 

1.360 

1.180 

1.028 

1.746              J 

..939 

1.344 

1.296 

1.626 

1.260 
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CASE      3      (C8NT) 
Yd/1)  YCI/2)  Y(I/3)  Y(I*4)  Y(Ii5>  Y(I/6) 


67 

1.203 

1.601 

1.525 

1.610 

1.575 

1.151 

68 

1.431 

1.050 

1.242 

1.270 

1.H6 

1.874 

69 

1.286 

1.170 

1.939 

1.011 

1.486 

1.768 

70 

1.C89 

1.497 

1.623 

1.972 

1.532 

1.876 
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CASE      3      (C©NT) 


X(J/1) 


X(J/2) 


XU,3) 


X(J/4) 


X(J,5) 


X(J/6) 


•  538 

►  236 

•  840 

.990 

.761 

.730 

•  063 

►  843 

.799 

►  273 

.804 

.882 

.251 

.784 

.924 

.409 

.471 

.011 

•  466 

.538 

.499 

.784 

.751 

.413 

•  293 

.315 

.990 

.536 

.259 

.084 

•  023 

.178      , 

.375 

.947 

.090 

.274 

•  239      , 

.141 

.505 

.949 

►  500 

.762 

•  402 

.277 

.659 

►  664 

►  139 

.676 

•  335      , 

.318       < 

.463       « 

.551 

.079 

.469 

•  573       , 

,189       , 

,935       « 

►  751 

►  155 

►  920 

•  971 

.905       , 

,657 

>175 

,275 

►  628 

.949       , 

,428       i 

,674       < 

,400 

,400 

►  171 

.967       , 

,284 

,305       i 

,450 

,688 

,055 

.928       , 

,429 

,543       < 

,477 

,709 

.642 

.112 

,149 

884 

,582 

,252 

.825 

•  694       , 

,080 

420 

»715 

,068 

,300 

•  171 

265 

371 

528 

462 

.282 

.543       , 

909 

612 

064 

572 

,150 

.233 

987 

346       < 

305 

748 

.877 

.439       , 

327       , 

231 

763 

198 

816 

.182       , 

594       , 

759 

751 

420 

027 

.334       . 

007 

058 

455 

079 

571 
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CASE      3       (C6NT) 
X(J/1)  X(J/2)  X(J,3)  X(J/4) 


X(J/5) 


X(J/6) 


#366 

.942       « 

.747 

►  010 

,371 

,237 

.067 

.288       , 

►  326      < 

,129 

,143 

,353 

•  304       « 

.483       < 

.219 

,195 

,713 

,821 

.000      < 

,257 

.693       < 

>116 

,194 

,278 

•  724       , 

,505 

.323       « 

,180 

,806 

,718 

.702 

,384       , 

,307 

,345 

,351 

,636 

•  506       « 

,967 

,298 

,029 

,102 

,656 

•  766       < 

,787       « 

,584 

,131 

,092 

,119 

•  038       « 

,807 

,826 

,044 

,988 

,259 

•  643       « 

,534       < 

,796 

,098 

,963 

,982 

.296       « 

,412 

,109 

,692 

,836 

,105 

.337       , 

,672      < 

,070 

,882 

,448 

,305 

•  398 

,659       , 

,515      « 

,768 

,378       , 

,627 

.476 

,091 

,688 

,552 

,897 

,989 

.477 

,892       < 

,514 

,633 

,048       i 

,520 

•  166       , 

,004 

,062 

,839 

,412 

,591 

•  596       , 

,112 

,665 

,498 

,376 

,041 

.366       , 

,208 

,867 

,493 

,592 

,738 

.060       « 

,130       < 

,542 

,230 

,975 

,796 

.433       , 

,373 

,658 

,486       < 

,854 

,466 

.955 

,541 

,412 

,861 

,513 

,525 

.112 

,364       , 

,799 

,807 

,380 

,097 
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CASE      3      (CONT) 
X(Jil)  X(J/2)  X(J,3)  X(J*4) 


X(J/5) 


X(J/6) 


45      t 

.325 

.321       « 

,120 

.010 

.338 

►  971 

46       < 

.134 

.692      « 

.961 

.127 

.770 

.188 

47       < 

.767 

.734       ( 

.132 

,436 

.145 

.476 

48       < 

.234       < 

.679 

.315 

►  683 

.080 

.980 

49       i 

.318 

.306       < 

.583       < 

,455       « 

.137 

.658 

50 

.890 

.423       « 

.538 

.126 

.138 

.626 

51 

.621       , 

.712 

.449       « 

.148       < 

.108 

.278 

52 

.421       , 

.105       , 

.245       * 

.284       « 

.934 

.728 

53      « 

.718 

,146 

.622 

.702 

.523 

.575 

54      i 

.826       , 

,635 

.266 

.808 

.058 

.599 

55 

.940       « 

.189       i 

.933       < 

.178 

.540       . 

.429 

56 

.636       ■ 

.978 

.629       < 

,380 

.292 

.874 

57       , 

.780       . 

,202      ( 

,703 

,072 

.150 

.824 

58       i 

.981       < 

,188       , 

,637 

.680       i 

.416 

.818 

59 

.884       ( 

,289 

,694       i 

•  340       < 

.876 

.109 

60 

.742 

,621       < 

.901 

•  104 

.638       < 

.532 

61 

.165 

,094       , 

,425 

.477 

.698 

.248 

62 

.558 

,528 

,688 

.278 

.953 

►  264 

63       . 

.570       , 

,175 

,717 

.916 

.320 

►  253 

64 

.511 

,066 

,864       < 

►  498       « 

►  872 

►  186 

65       , 

.559       , 

,044       . 

,175 

•  705 

►  016 

►  562 

66       , 

.403       < 

,931 

,299 

►  654       < 

.600 

►  733 
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CASE      3      (C6NT) 
J  X(J/1)  X(J/2)  X(J,3)  X(J/4)  X(J,5>  X(Ji6) 


67 

.463 

•  481 

•  116 

•  074 

.870 

.213 

68 

.739 

•  094 

.651 

•  967 

.349 

.173 

69 

•  504 

.663 

•  638 

•  712 

.513 

.822 

70 

.212 

•  056 

•  592 

.658 

.391 

.566 
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..*.«v  « 


CASE      4      DATA 


Yd/1) 


Y(I*2> 


Y(I#3) 


Y(I*4) 


1 

.2500 

-1.2000 

.7000 

-1. 

.0000 

2     -1. 

.2000 

-.0500 

-1. 

.2000 

^ 

.0500 

3 

.5000 

-1.0000 

-  < 

.5000 

-1. 

.0000 

4      1. 

.0000 

-.7000 

-< 

.6000 

1< 

.0000 

5      1. 

.0000 

.5000 

-1. 

.0000 

.8000 

6      It 

.2000 

.7000 

•  i 

.7000 

-1. 

►  1000 

7 

.5000 

1.0000 

.6000 

1. 

.2000 

8 

.5250 

1.0540 

1. 

.1500 

.0550 

9       1< 

.1000 

1.1000 

.2000 

.2000 

10       1< 

.1000 

-1.1000 

•  i 

.2000 

.2000 

11      -1< 

.1000 

1.1000 

.2000 

m   ^ 

.2000 

12      -1< 

.1000 

-1.1000 

m   , 

.2000 

«■  i 

.2000 

13       1, 

.5000 

.5000 

.8000 

.6000 

m          i. 

.5000 

-.5000 

< 

.7500 

m   < 

.7500 

15 

.0470 

•1.5000 

—  < 

.3000 

.9400 

16      -1< 

.2000 

.5000 

1. 

.0000 

.5000 

17     »li 

>2000 

-.5000 

1. 

.0000 

-< 

.5000 

18       1, 

0000 

-.4000 

-1« 

4000 

.4000 

19      1, 

.0000 

-.4000 

-1« 

4000 

m  < 

.4000 

20       i 

2000 

-.1000 

1. 

.0000 

1. 

.3000 

21 

.6820 

-.3590 

•  i 

.1320 

-1< 

.8380 

22 

7650 

1.4370 

.6250 

-8080 
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CASE      4       <C6NT) 
Y(N1)  Y(Ii2) 


Y(Ii3) 


Y(I/4) 


23 

►  5220 

-.4680       -i 

►  2980 

-1. 

.8480 

24 

•  4 

•  7580 

--.5810       1« 

.7140 

►  3730 

25 

«1< 

►  9940 

.0860 

►  1110 

m   ( 

►  0550 

26 

►  7020 

••7940       1. 

.6710 

►  2390 

27 

►  3690 

-.8250       1. 

.6220 

►  5890 

28 

V  ( 

.8440 

-.2490 

.0820 

-1. 

►  7860 

29 

"  < 

►  7320 

.0710 

.4260 

-1< 

.7740 

30 

«•  | 

►  7350 

.0330 

.0700 

-1. 

►  8580 

31 

-1, 

►  9130 

.3420 

.1820 

.2230 

32 

m  i 

.8180 

1.4470       -< 

►  9770 

.0940 

33 

m  i 

.3800 

1.4660      -1, 

.1350 

-, 

.6260 

34 

•  < 

►  5510 

1.5800 

.9870 

.4580 

35 

«1« 

.9410 

-.3660        , 

►  1630 

.0220 

36 

1130 

1.5680      -1. 

.0680 

.3450 

37 

0980 

1.5600 

.3780 

-1. 

.1680 

38 

"1" 

9060 

.2180 

.4960 

.1410 

39 

-1« 

.0050 

-.4050       1< 

.6320 

►  1690 

40 

-1« 

9770 

-.2900 

.0450 

•  i 

.0510 
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CASE 


(CQNT) 


X(Jil) 


X(J/2) 


XU/3) 


X(J/4) 


1 

►  1370 

►  6680 

,5700 

,4500 

2 

►  2450 

►  7640 

►  5690 

,0990 

3 

►  1140 

►8740        < 

►  2970 

►  3620 

4        « 

►  4360 

►5580        i 

.4460        < 

►  2870 

5 

►  1970 

►  5910 

,4740 

,5380 

6 

►  0770 

►  9320 

.2750 

►  1710 

7 

►  2170 

►  6450 

.6260 

►  1710 

8        < 

.0970       1 

►  5710 

►7950        < 

►  1120 

9 

►  2930 

►6520        « 

►  2490 

►  2060 

10 

►  2080 

►  7060 

►  4710 

►  3580 

11 

rl990 

►8330        « 

►  3580 

►  1330 

12        , 

►  2670 

►7540       i 

,3010 

►  2510 

13 

►0280      -« 

►  9320 

►  3580 

►  0480 

14 

.4860       -< 

.7850 

.1990 

►  2100 

15 

.2540 

6090        . 

►2680       -< 

.5980 

16        , 

1630 

.8880        < 

►2430        < 

►  2220 

17 

.4100        < 

.3490        , 

.3970 

.3460 

18 

2020       «•< 

.7230 

►  5960 

►  1500 

19 

.3550        i 

.5350 

.6690 

.3320 

20 

.2280       < 

5960        , 

.4750       -, 

.5370 
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CASE      5      DATA 


Ydil) 


Y(J*2> 


YII/3) 


Y(I*4) 


1 

5« 

.1000 

3< 

►  5000 

1.4000 

.2000 

2 

4« 

.9000 

3< 

►  0000 

1.4000 

.2000 

3 

4< 

►  7000 

3. 

►  2000 

1.3000        « 

►  2000 

4 

4< 

►  6000 

3. 

.1000 

1.5000        « 

►  2000 

5 

5< 

►  0000 

3- 

►  6000 

1.4000        « 

.2000 

6 

5< 

►  4000 

3« 

►  9000 

1.7000        « 

.4000 

7 

<H 

►  6000 

3< 

►  4000 

1.4000        < 

.3000 

8 

5< 

►  0000 

3. 

►  4000 

1.5000        , 

.2000 

9 

4< 

•  4000 

2. 

►  9000 

1.4000        . 

.2000 

10 

4< 

►  9000 

3< 

►  1000 

1.5000 

,1000 

11 

5< 

►  4000 

3< 

.7000 

1.5000 

.2000 

12 

4< 

►  8000 

3. 

►  400C 

1.6000 

.2000 

13 

4< 

►  8000 

3< 

►  0000 

1.4000        < 

.1000 

H 

4< 

►  3000 

3i 

►  0000 

1.1000 

.1000 

15 

5« 

►  8000 

4< 

►  0000 

1.2000 

.2000 

16 

5« 

►  7000 

4< 

►  4000 

1.5000        < 

►  4000 

17 

5. 

.4000 

3< 

►  9000 

1.3000 

.4000 

18 

5« 

►  1000 

3. 

►  5000 

1.4000        . 

.3000 

19 

5. 

►  7000 

3« 

►  8000 

1.7000        < 

.3000 

20 

5< 

►  1000 

3. 

►  8000 

1.5000 

.3000 

21 

5« 

►  4000 

3< 

►  4000 

1.7000        « 

.2000 

22 

5 

►  1000 

3. 

►  7000 

1.5000 

.4000 
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CASE      5      (C9NT) 
Y(Iil)  Y(li2)  Y(I/3) 


Y(I*4) 


23 

4« 

.6000 

3< 

.6000 

1. 

.0000 

.2000 

24 

5. 

.1000 

3. 

.3000 

1. 

.7000 

.5000 

25 

4« 

.8000 

3< 

.4000 

1. 

.9000 

.2000 

26 

5< 

.9000 

3 

.2000 

4 

►  8000 

1.8000 

27 

6< 

.1000 

2. 

.8000 

4. 

.0000 

1»3000 

28 

6< 

►  3000 

2« 

.5000 

4. 

•  9000 

1.5000 

29 

6< 

.1000 

2< 

.8000 

4. 

.7000 

1.2000 

30 

6< 

.4000 

2. 

.9000 

4. 

.3000 

1.3000 

31 

7« 

.0000 

3< 

.2000 

4. 

.7000 

1.4000 

32 

6« 

.4000 

3< 

.2000 

4< 

.5000 

1.5000 

33 

6. 

•  9000 

3« 

.1000 

4. 

.9000 

1.5000 

34 

5. 

5000 

2. 

•  3000 

4. 

.0000 

1.3000 

35 

6« 

5000 

2< 

.8000 

4, 

.6000 

1.5000 

36 

5« 

7000 

2. 

.8000 

4. 

.5000 

1.3000 

37 

6« 

3000 

3« 

3000 

4, 

.7000 

1.6000 

38 

4. 

9000 

2« 

.4000 

3. 

.3000 

1.0000 

39 

6« 

-6000 

2< 

.9000 

4. 

.6000 

1.3000 

40 

5< 

2000 

2< 

.7000 

3. 

.9000 

1.4000 

41 

5< 

0000 

2« 

0000 

3< 

.5000 

1.0000 

42 

5. 

9000 

3« 

0000 

4. 

.2000 

1.5000 

43 

6. 

0000 

2. 

2000 

4, 

.0000 

1.0000 

44 

6. 

1000 

2< 

9000 

4« 

7000 

1.4000 
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CASE      5 

(C6NT) 

I 

YCI/1) 

Yd/2) 

YU#3> 

YU/4) 

45 

5.6000 

2.9000 

3.6000 

1.3000 

46 

6.7000 

3.1000 

<f.4000 

1.4000 

47 

5.6000 

3.0000 

4.5000 

1.5000 

48 

5.8000 

2.7000 

4.1000 

1.0000 

49 

6.2000 

2.2000 

4.5000 

1.5000 

50 

5.6000 

2.5000 

3.9000 

1.1000 

• 
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CASE      5       (CSNT) 


XiJtl) 


X(J/2> 


X(J,3) 


X(J/4) 


1 

6< 

.5000 

3< 

.2000 

5. 

1000 

2.0000 

2 

6< 

.4000 

2« 

.7000 

5* 

3000 

1»9000 

3 

6. 

»8000 

3< 

.0000 

5. 

.5000 

2.1000 

4 

5< 

.7000 

2< 

.5000 

5. 

.0000 

2.0000 

5 

5< 

.8000 

2. 

.8000 

5. 

.1000 

2.4000 

6 

61 

.4000 

3« 

.2000 

5. 

.3000 

2.3000 

7 

6< 

.5000 

3. 

.0000 

5. 

.5000 

1.8000 

8 

7« 

.7000 

3< 

.8000 

6< 

.7000 

2.2000 

9 

7< 

.7000 

2< 

.6000 

6< 

.9000 

2.3000 

10 

6« 

.0000 

2< 

.2000 

5. 

.0000 

1.5000 

11 

6< 

.3000 

3< 

.3000 

6, 

.0000 

2.5000 

12 

5< 

.8000 

2« 

.7000 

5. 

.1000 

1.9000 

13 

7« 

•  1000 

3. 

.0000 

5, 

.9000 

2.1000 

14 

6< 

.3000 

2< 

.9000 

5, 

.6000 

1.8000 

15 

6< 

.5000 

3< 

.0000 

5, 

.8000 

2.2000 

16 

7. 

.6000 

3« 

.0000 

6, 

.6000 

2.1000 

17 

<H 

►  9000 

2« 

.5000 

4. 

.5000 

1.7000 

18 

7« 

.3000 

2< 

.9000 

6« 

.3000 

1.8000 

19 

6i 

.7000 

2- 

.5000 

5. 

.8000 

1.8000 

20 

7 

.2000 

3< 

.6000 

6< 

.1000 

2.5000 

21 

6 

.9000 

3< 

.2000 

5« 

.7000 

2.3000 

22 

5 

.6000 

2< 

.8000 

4 

r9000 

2.0000 
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CASE   5   (C6NT) 
X(J#1)       X(J*2> 


X(J/3) 


X(J,4) 


23 
2* 
25 


7.7000  2.8000 
6.3000  2.7000 
6.7000      3.3000 


6.7000 
*f.9000 

5.7000 


2.0000 
1.8000 
2.1000 
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CASE      6      DATA 


I         YU,1>    Yd*2)    Y(I/3)    Y(Ii4)    Yd, 5)    Yd/6)    Yd/7)    Yd/8) 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

15 
16 
17 
18 
19 
20 
21 
22 


160. 

►     5. 

.1 

10. 

.0 

28. 

155. 

5. 

.2 

8. 

.0 

27. 

130. 

►     5. 

.0 

10. 

.0 

26. 

161. 

►     5. 

.0 

10. 

.0 

27. 

135. 

5. 

.0 

12« 

.0 

27. 

165. 

>     5< 

.0 

11« 

.0 

23. 

150. 

.     5. 

•  0 

9. 

.0 

29. 

148. 

4. 

.8 

8. 

.0 

26. 

150. 

.     4. 

.5 

7. 

.0 

31. 

120. 

4. 

.0 

6. 

.0 

33. 

120. 

5. 

.1 

8. 

.0 

32. 

190. 

>     4. 

»2 

8. 

.0 

30. 

100. 

•     4< 

►  4 

9. 

.0 

35. 

150. 

►     4. 

►  0 

7. 

►  0 

29. 

90. 

►     4< 

>6 

8. 

►  0 

30. 

120. 

.     4. 

.7 

7. 

.0 

35. 

200. 

►     4< 

.3 

9< 

.0 

30. 

120. 

>     4. 

►  1 

8< 

.0 

28. 

105. 

5 

.0 

7< 

►  0 

27. 

210- 

►     5 

.2 

8 

.0 

26. 

90 

►     4 

►  0 

10 

.0 

25. 

110 

»     5 

.2 

11 

.0 

25. 

0 

70. 

450. 

4. 

.10 

21. 

0 

0 

85. 

400. 

3. 

.12 

17. 

0 

0 

75. 

560. 

►    2. 

»96 

19. 

.0 

0 

70. 

665. 

>    4. 

.04 

21. 

0 

0 

88. 

570. 

.    3. 

.05 

23. 

0 

0 

95. 

675. 

3. 

.72 

24. 

•  0 

0 

90. 

580. 

.    3. 

.83 

25. 

.0 

0 

85. 

390. 

3. 

.75 

28. 

.0 

0 

60. 

435. 

,    3. 

.05 

20« 

.0 

0 

55. 

440. 

2. 

.78 

18. 

.0 

0 

56. 

650. 

2. 

.89 

23. 

-0 

0 

55. 

640. 

3. 

.46 

27. 

.0 

0 

48. 

430. 

2 

.98 

26. 

.0 

0 

65. 

650. 

2. 

►  84 

20. 

.0 

0 

70. 

655. 

2« 

.96 

17. 

.0 

0 

67. 

645. 

2. 

.78 

20. 

.0 

0 

62. 

660. 

3. 

.20 

23. 

.0 

0 

63. 

530. 

3 

►  05 

19. 

.0 

0 

64. 

435. 

3 

►  55 

27 

.0 

0 

67. 

440. 

>    3 

►  10 

20 

.0 

0 

68. 

430 

3 

►  41 

21 

.0 

0 

60. 

530 

►    2 

►  97 

25. 

.0 

82 


CASE      6       (C9NT) 
I         Ylhl)    Yd, 2)    YU>3)    Yd/4)    Yd*5)    Yd, 6)    Yd, 7)    Y(Ii8) 


23 
24 
25 
26 

27 

28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 


100. 

.     4.3 

9.0 

25. 

70. 

»     4.5 

8.0 

23. 

100. 

.     4*8 

9.0 

27. 

130. 

»     5.2 

9.0 

25. 

90. 

,           4.5 

11.0 

37. 

80 

►    .  4.6 

10.0 

32. 

95. 

4.9 

10.0 

25. 

70. 

>     4.4 

12.0 

30. 

80. 

»     4.2 

13.0 

36. 

95. 

5.1 

15.0 

31. 

100. 

i     4.6 

11.0 

24. 

95. 

4.8 

10.0 

27. 

85. 

4.7 

12.0 

25. 

70. 

4.8 

11.0 

26. 

80. 

•     5.4 

10.0 

21. 

85. 

5.5 

13.0 

33. 

.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 
.0 


70. 
64, 
65. 
70. 
74, 
78. 
82. 
85. 
69. 
70. 
76. 
74. 
73. 
78. 
80. 
81. 


440 
450 
355 
380 
350 
450 
260 
262 
260 
265 
270 
355 
360 
365 
370 
355 


62 

17. 

.0 

40 

19. 

.0 

10 

24. 

.0 

80 

29. 

.0 

95 

31. 

.0 

12 

32 

.0 

11 

23. 

>0 

22 

24. 

►  0 

35 

26. 

.0 

40 

29. 

.0 

54 

30. 

.0 

63 

31. 

.0 

17 

29. 

.0 

18 

30. 

.0 

90 

32. 

.0 

12 

35. 

•  0 

83 


CASE      6       (C6NT) 


X(J,1)    X(J#2)    X(J/3>    X(Ji4)    X(Ji5>    X(Ji6)    X(J,7)    X(J#8> 


1 
2 
3 

4 
5 
6 
7 

8 
9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 


200. 

.     3. 

>4 

10. 

>0 

24. 

260. 

.     3. 

>1 

8. 

»0 

21. 

195< 

3. 

.0 

9. 

.0 

20. 

195. 

3. 

.2 

9. 

.0 

19. 

220. 

»     3. 

.3 

10. 

.0 

24. 

220. 

►     3. 

.0 

8. 

►  0 

25. 

190. 

3< 

.4 

9. 

»0 

26. 

285. 

3. 

»0 

11. 

.0 

19. 

300. 

i     3. 

»0 

9. 

>0 

20. 

225. 

3. 

•  0 

10. 

.0 

22. 

260. 

.     3. 

>4 

8. 

.0 

22. 

280. 

.     3- 

•  0 

8. 

.0 

20. 

300. 

3< 

•  4 

10. 

.0 

20. 

310. 

3. 

.0 

11. 

►  0 

19. 

290. 

•     3. 

»1 

12. 

.0 

26. 

260. 

.     3. 

.0 

8. 

►  0 

22. 

290. 

►     3. 

.3 

9. 

.0 

25. 

160. 

»     3. 

.1 

11. 

.0 

20. 

240. 

►     3. 

.5 

11. 

>0 

20. 

195. 

►     3 

►  1 

8. 

.0 

21. 

290. 

►     3- 

.4 

10. 

►  0 

19. 

210 

.     3 

►  5 

9- 

.0 

22. 

98. 
110. 
105. 
110. 

95. 

90. 

96. 
100. 
1C2. 
lC5t 

97. 
112. 
108. 
106. 

94. 

98. 
100. 

79. 

88. 

81. 

94. 

96. 


210 
220 
130 
010 
205 
210 
070 
990 
120 
985 
090 
200 
835 
055 
240 
015 
010 
170 
990 
975 
860 
950 


10 

40. 

0 

89 

36. 

0 

98 

44. 

.0 

10 

36. 

0 

54. 

29. 

.0 

78 

34. 

.0 

36 

30. 

-0 

23 

31. 

.0 

75 

40. 

.0 

80 

25. 

.0 

94 

38< 

.0 

70 

41. 

.0 

20 

45. 

.0 

45 

37< 

>0 

50 

41. 

►  0 

60 

29. 

.0 

56 

32. 

.0 

26 

34. 

.0 

40 

42 

.0 

58 

44. 

-0 

40 

40. 

.0 

70 

47. 

»0 

84 


CASE      6       (CeNT) 
J         X(Jil>    X(J/2>    X(J/3)    X(J/4>    X(J/5)    X(J,6>    XU/7)    X(J/8) 


23 

180 

24 

205 

25 

215 

26 

270 

27 

290 

28 

320 

29 

210 

30 

210 

31 

185 

32 

200 

33 

170 

34 

140 

35 

110 

36 

100 

37 

95 

38 

65 

39 

55 

40 

70 

4i 

85 

42 

115 

43 

95 

44 

120 

•     3- 

.0 

11« 

.0 

22. 

.0 

2< 

.9 

11< 

.0 

23. 

.0 

•     3< 

»4 

8< 

.0 

21. 

0 

3. 

.1 

8< 

.0 

20- 

.0 

3. 

.0 

9< 

»0 

23. 

.0 

3. 

.2 

10< 

.0 

19. 

.0 

3< 

.0 

9< 

.0 

18. 

.0 

3. 

.0 

8. 

.0 

21. 

»0 

3< 

►  4 

9. 

.0 

25. 

.0 

3. 

.2 

8< 

.0 

26< 

>0 

.     2. 

►  9 

9< 

►  0 

20. 

rO 

3. 

.0 

9< 

.0 

20. 

.0 

•     4< 

.9 

9< 

.0 

22. 

'0 

5« 

.6 

8< 

»0 

19. 

.0 

4. 

,9 

8. 

.0 

24. 

•  0 

6. 

.2 

9. 

.0 

30< 

iO 

5. 

.0 

10« 

.0 

27. 

.0 

5. 

.3 

7< 

►  0 

28. 

»0 

<H 

.9 

11- 

.0 

19. 

.0 

5. 

>0 

10« 

►  0 

21. 

.0 

4. 

.8 

8< 

.0 

27 

►  0 

.     6 

.1 

9 

.0 

24 

»0 

97. 

90. 
100. 
111. 
102. 

87. 
112. 

95» 

96. 

98. 

95. 

98. 
130. 
128. 
124. 
134. 
128. 
118. 
117. 
122. 
114. 
120. 


990 

805 
700 
170 
350 
160 
010 
190 
055 
980 
095 
990 
220 
216 
218 
200 
205 
204 
206 
198 
228 
244 


83 

28< 

-0 

25 

38. 

.0 

05 

32. 

0 

60 

31. 

>0 

42 

33. 

•  0 

64 

38. 

•  0 

81 

45. 

,0 

00 

40. 

.0 

18 

41. 

►  0 

46 

38- 

.0 

15 

40. 

-0 

25 

32. 

.0 

51 

19. 

.0 

60 

14 

.0 

45 

10 

•  0 

61 

12 

.0 

62 

11 

.0 

35 

12 

.0 

90 

14 

►  0 

25 

17 

.0 

65 

11 

.0 

82 

18 

.0 

85 


CASE   7   DATA 


I    Yl   Y2    Y3    Y4    Y5   Y6   Y7   Y8   Y9   YlO  Yll  Yl2   Y13 


1    i 

,   365 

152 

90 

4 

2 

4 

3 

1 

1 

5 

250 

39 

2    1 

I   517 

126 

70 

2 

3 

2 

2 

2 

2 

326 

38 

3    1 

I   343 

180 

86 

1 

3 

2 

2 

3 

2 

278 

42 

*»    3 

L   432 

110 

65 

1 

3 

2 

4 

4 

2 

278 

42 

5    1 

I   376 

115 

80 

1 

3 

2 

1 

3 

2 

278 

41 

6    1 

L   410 

110 

70 

1 

3 

2 

2 

3 

2 

255 

43 

7    j 

L   408 

105 

65 

1 

3 

2 

2 

3 

2 

278 

50 

8    1 

L   410 

150 

100 

2 

4 

2 

2 

4 

2 

278 

44 

9    1 

I   315 

118 

76 

2 

3 

2 

2 

4 

5 

315 

28 

10    1 

L   546 

150 

70 

1 

3 

2 

2 

3 

1 

212 

47 

11    1 

1   437 

120 

70 

2 

3 

2 

2 

4 

? 

372 

39 

12   j 

I   506 

106 

72 

2 

3 

2 

2 

3 

145 

37 

13    J 

I   338 

110 

70 

2 

3 

2 

1 

3 

278 

39 

H   : 

I   446 

120 

70 

1 

2 

2 

1 

4 

278 

50 

15      : 

1   426 

135 

80 

1 

3 

2 

1 

4 

278 

48 

16      : 

I   410 

140 

85 

1 

3 

2 

2 

4 

250 

35 

17 

L   483 

120 

70 

2 

3 

2 

2 

4 

278 

48 

is   : 

3   377 

114 

70 

2 

2 

2 

2 

2 

4 

5 

350 

38 

19   : 

3   410 

130 

85 

2 

1 

3 

2 

2 

4 

2 

264 

39 

20   : 

3   414 

120 

80 

4 

2 

1 

2 

1 

3 

5 

246 

49 

21 

I   369 

140 

90 

4 

1 

3 

2 

1 

4 

1 

278 

38 

22 

i   335 

132 

82 

3 

1 

3 

1 

1 

4 

1 

278 

38 

86 


I    Yl   Y2    Y3 


CASE   7   (C8NT) 
Y4    Y5   Y6   Y7   Y8   Y9   YlO  Yll  Y12   Y13 


23    3 

,   410 

140 

80    2 

J    1 

3 

1 

2 

4 

2 

268 

40 

24    J 

I   422 

152 

92    3 

2 

2 

2 

2 

3 

2 

278 

49 

25    1 

I   517 

110 

80    3 

2 

3 

2 

2 

4 

5 

250 

53 

26    1 

.   250 

130 

70    3 

2 

3 

2 

2 

4 

2 

315 

45 

27    3 

.   466 

114 

74    3 

2 

2 

2 

4 

3 

250 

38 

28    3 

1   366 

120 

90    3 

2 

2 

2 

4 

3 

278 

47 

29    3 

I   271 

110 

70    3 

3 

2 

2 

4 

5 

295 

36 

30    J 

I   467 

130 

70    « 

3 

2 

2 

4 

5 

278 

41 

31    1 

l   372 

128 

92    < 

3 

2 

2 

4 

5 

353 

43 

32    3 

L   300 

140 

75    l 

\        2 

3 

2 

2 

4 

5 

200 

38 

33    J 

I   348 

110 

80    « 

2 

1 

2 

4 

2 

335 

40 

34    3 

I   288 

130 

80    1 

4 

2 

2 

4 

2 

264 

47 

35    3 

I   477 

115 

80    3 

2 

2 

2 

k 

2 

245 

43 

36    3 

L   406 

140 

110     t 

J    2 

2 

2 

2 

4  . 

5 

300 

33 

37    3 

I   489 

110 

64    1 

I    2 

2 

2 

2 

3 

5 

278 

42 

38    3 

L   406 

130 

80    ' 

f    2 

4 

2 

2 

3 

5 

278 

51 

39    3 

L   367 

135 

85    3 

I    2 

3 

2 

1 

3 

3 

350 

40 

40      : 

I   444 

125 

80    3 

I    2 

3 

2 

1 

3 

1 

302 

45 

4i      : 

I   402 

122 

70   : 

I    1 

1 

2 

2 

3 

1 

278 

^6 

42      : 

I   462 

120 

78   : 

L    2 

3 

1 

2 

4 

1 

278 

41 

43      : 

I   511 

140 

95    l 

♦    1 

2 

2 

1 

3 

1 

350 

45 

44      : 

I   410 

140 

90    J 

*    1 

2 

2 

1 

3 

? 

245 

52 

87 


CASE   7   (C9NT) 
I    Yl   Y2    Y3    Y4    Y5   Y6   Y7   Y8   Y9   Y10  YU  Y12   Y13 


45 

1 

446 

90 

60 

1    1 

4 

1 

1 

3 

2 

202 

40 

46 

1 

468 

110 

70 

1    2 

3 

2 

1 

4 

1 

314 

37 

47 

1 

410 

130 

80 

1    2 

3 

2 

2 

4 

5 

351 

40 

48 

1 

410 

140 

88 

1    2 

3 

2 

2 

4 

1 

242 

42 

49 

I 

355 

120 

70 

1    1 

3 

2 

2 

4 

3 

366 

23 

50 

2 

394 

116 

80 

1    2 

1 

2 

2 

4 

2 

283 

42 

51 

2 

312 

100 

60 

1    2 

1 

2 

2 

4 

2 

310 

47 

52 

1 

383 

100 

60 

1    2 

1 

2 

2 

4 

2 

310 

47 

53 

1 

489 

140 

80 

1    2 

2 

2 

1 

4 

2 

278 

46 

54 

1 

410 

130 

82 

1    2 

2 

2 

1 

4 

2 

250 

41 

55 

1 

404 

140 

98 

1    2 

2 

2 

I 

4 

2 

278 

43 

56 

1 

378 

140 

90 

1    2 

3 

1 

2 

3 

1 

267 

43 

57 

1 

369 

130 

70 

2    1 

4 

2 

1 

3 

1 

278 

44 

58 

1 

388 

133 

90 

3    1 

3 

2 

3 

4 

5 

288 

43 

59 

1 

483 

110 

80 

4    1 

2 

2 

2 

4 

2 

277 

39 

60 

3 

461 

140 

90 

1    2 

4 

1 

1 

2 

1 

330 

52 

88 


CASE 


(C9NT) 


XI   X2    X3 


X4 


X5   X6   X7   X8   X9   X10  Xll  Xl2   X13 


1 

1 

404 

160 

78 

4 

2 

1 

2    1 

4    1 

200 

67 

2 

1 

418 

118 

76 

1 

1 

3 

3    I 

1    2 

298 

43 

3 

1 

288 

128 

86 

1 

2 

4 

3    1 

1    2 

209 

37 

4 

2 

349 

170 

110 

3 

2 

1 

1    4 

1    5 

343 

40 

5 

1 

396 

118 

72 

3 

1 

2 

3    1 

1    2 

221 

37 

6 

3 

411 

180 

120 

1 

2 

4 

1    1 

1    2 

196 

52 

7 

1 

443 

160 

90 

4 

2 

3 

3    1 

1    2 

271 

37 

8 

1 

479 

140 

90 

1 

2 

2 

1    1 

1    1 

260 

40 

9 

3 

388 

140 

90 

1 

2 

3 

1    1 

1    5 

281 

40 

10 

2 

445 

130 

88 

4 

2 

4 

3    1 

1    1 

215 

42 

11 

1 

389 

110 

78 

3 

2 

4 

3    I 

1    2 

249 

32 

12 

1 

406 

114 

70 

2 

2 

3 

1    1 

1    5 

248 

40 

13 

1 

431 

100 

60 

4 

2 

3 

1    1 

1    1 

224 

34 

14 

2 

413 

144 

86 

4 

2 

4 

3    1 

1    1 

247 

41 

15 

1 

411 

no 

50 

4 

1 

4 

1    I 

1    2 

240 

40 

16 

1 

406 

110 

76 

1 

1 

3 

3    1 

4    2 

248 

23 

17 

1 

444 

120 

80 

1 

2 

1 

3    1 

1    1 

151 

42 

18 

2 

375 

120 

70 

1 

2 

3 

3    1 

1    3 

248 

41 

19 

1 

383 

120 

80 

2 

1 

3 

3    1 

1    1 

256 

36 

20 

2 

439 

114 

80 

1 

2 

4 

3    1 

1    2 

2^8 

40 

21 

2 

489 

104 

70 

1 

2 

3 

3    1 

3    1 

221 

44 

22 

1 

387 

118 

90 

1 

1 

4 

3    1 

1    2 

279 

43 

89 


XI   X2    X3 


CASE   7   (CSNT) 
X4    X5   X6   X7   X8   X9   XlO  Xll  X12   X13 


23 

3 

351 

120 

90 

4 

2 

3 

3   2 

5 

5 

258 

36 

24 

1 

350 

110 

66 

1 

1 

3 

3    1 

1 

3 

248 

41 

25 

3 

536 

124 

84 

4 

2 

2 

3    1 

1 

5 

282 

33 

26 

2 

432 

110 

70 

1 

1 

4 

1    1 

2 

1 

240 

55 

27 

3 

393 

150 

90 

3 

1 

3 

1    1 

1 

5 

245 

46 

28 

1 

369 

132 

80 

1 

2 

2 

3    1 

1 

1 

237 

40 

29 

2 

418 

120 

80 

3 

2 

3 

3    1 

3 

2 

236 

40 

30 

383 

126 

84 

1 

1 

4 

3    1 

0 

1 

221 

36 

31 

325 

135 

80 

1 

2 

4 

3    1 

5 

5 

175 

30 

32 

395 

120 

78 

4 

1 

4 

1    1 

1 

178 

40 

33 

375 

117 

75 

1 

1 

4 

3    1 

1 

229 

39 

34 

432 

142 

80 

4 

2 

4 

3    1 

1 

235 

43 

35 

267 

120 

84 

4 

2 

4 

3    1 

1 

213 

40 

36 

276 

158 

100 

1 

2 

4 

3    1 

2 

167 

41 

37 

489 

188 

70 

1 

2 

3 

3    1 

3 

248 

44 

38 

369 

120 

84 

1 

1 

4 

3    1 

1 

294 

43 

39 

398 

150 

90 

4 

2 

4 

3    1 

1 

223 

42 

t 

40 

442 

118 

70 

3 

2 

2 

3    1 

5 

240 

44 

41 

394 

130 

80 

1 

2 

3 

3    1 

5 

288 

31 

42 

444 

120 

78 

3 

2 

2 

3    1 

2 

1 

256 

42 

43 

389 

140 

95 

1 

1 

k 

3    2 

2 

5 

273 

44 

44 

359 

130 

100 

1 

1 

2 

3    1 

1 

5 

195 

35 

90 


CASE   7   (C8NT) 
J    XI   X2    X3    X4    X5   X6   X7   X8   X9   XlO  Xll  X12   X13 


45 

2 

232 
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X       1 

1    5 
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34 

47 
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1    1 

1    1 
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39 
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1    1 

1    1 
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54 
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364 
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4 
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3    1 

1    5 
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33 
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140 

90 
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3    1 

1    5 

248 

36 
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APPENDIX  IV.   TABLE  OF  A  AND  B  VECTORS 

The  following  table  is  the  values  of  the  various  components 
of  the  A  and  B  vectors  used  in  the  seven  examples.   Only  six 
of  the  thirteeen  components  have  been  listed.   To  use  the 
table,  repeat  the  pattern  of  the  values  until  the  number  of 
components  equals  the  dimension  of  the  problem.   For  the 
remaining  components,  use  a  zero  value  until  one  has  filled 
the  thirteen  components. 


TABLE  OF  A  AND  B  VALUES 

DIMENSION 

COMPONENTS 

1 

2 

3 

4 

5 

6 

3 

A 

.50000 

.70711 

.50000 

.00000 

.00000 

.00000 

B 

.50000 

-.70711 

.50000 

.00000 

.00000 

.00000 

4 

A 

.50000 

.50000 

.50000 

.50000 

.00000 

.00000 

B 

.50000 

-.50000 

.50000 

-.50000 

.00000 

.00000 

5 

A 

.40825 

.50000 

.40825 

.50000 

.40825 

.00000 

B 

.40825 

-.50000 

.40825 

-.50000 

.40825 

.00000 

6 

A 

.40825 

.40825 

.40825 

.40825 

.40825 

.40825 

B 

.40825 

-.40825 

.40825 

-.40825 

.40825 

-.40825 

7 

A 

.35355 

.40825 

.35355 

.40825 

.35355 

.40825 

B 

.35355 

-.40825 

.35355 

-.40825 

.35355 

-.40825 

8 

A 

.35355 

.35355 

.35355 

.35355 

.35355 

.35355 

B 

.35355 

-.35355 

.35355 

-.35355 

.35355 

-.35355 

9 

A 

.31628 

.35355 

.31628 

.35355 

.31628 

.35355 

B 

.31628 

-.35355 

.31628 

-.35355 

.31628 

-.35355 

10 

A 

.31628 

.31628 

.31628 

.31628 

.31628 

.31628 

B 

.31628 

-.31628 

.31628 

-.31628 

.31628 

- .31628 

11 

A 

.28867 

.31628 

.28867 

.31628 

.28867 

.31628 

B 

.28867 

-.31628 

.28867 

-.31628 

.28867 

-.31628 

12 

A 

.28867 

.28867 

.28867 

.28867 

.28867 

.28867 

B 

.28867 

-.28867 

.28867 

-.28867 

.28867 

-.28867 

13 

A 

.26726 

.28867 

.26726 

.28867 

.26726 

.28867 

B 

.26726 

-.28867 

.26726 
92 

-.28867 

.26726 

-.28867 

APPENDIX  V.   TWO  DIMENSIONAL  PLOTS 

The  following  figures  represent  the  plots  of  the  computer 
printout.   The  values  of  the  A  and  B  vectors  are  printed  below 
the  plot  if  a  discriminant  function  was  determined.   The 
components  for  the  initial  projections  can  be  extracted  from 
Appendix  IV.   A  polaroid  print  is  also  included  for  each  case. 
This  allows  for  a  comparison  between  the  AGT  display  and  the 
.computer  plot  to  be  obtained. 
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Figure    IV-1-1  ,    CASE    1  ,    STEP    1 
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Figure    IV-1-2.      CASE    1  ,    STEP    2 
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Figure    IV-1-3.      Case    1 ,    STEP    3 
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Figure  IV-1-4.   CASE  1 ,  STEP  4 
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Figure    IV-1-5.      CASE    1 ,    STEP    5 
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Figure    IV-1-5P.      CASE    1,    STEP    5 
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Figure    IV-2-1 .      CASE    2,    STEP    1 
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Figure  IV-2-2 
A(l)  =-.0418 
A(2)  =  .9983 
A(3)  =-.0418 


CASE  2,  STEP  2 

B(l)  =  .7059 

B(2)  =  .0591 

B(3)  =  .7059 
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Figure  IV-2-3 
A(l)  =-.3240 
A(2)  =-.0005 
A(3)    =-.9461 


CASE  2,  STEP  3 
B(l)  =-.0248 
B(2)  =  .9997 
B(3)  =  .0079 
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Figure  IV-2-4. 
A(l)  =  .0081 
A(2)  =-.9998 
A(3)  =-.0161 


CASE  2,  STEP  5 
B(l)  =-.8980 
B(2)  =-.0002 
B(3)  =-.4400 
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Figure    IV-2-5. 
A(l)    =-.9122 
A(2)    =    .2671 
A(3)    =-.3107 


CASE  2,  STEP  6 

B(l)  =  .0887 

B(2)  =  .8690 

B(3)  =  .4867 
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Figure    IV-2-5P,    CASE    2,    STEP    6 
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Figure    IV-3-1  .      CASE    3,    STEP    1 
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Figure    IV-3-2.      CASE    3,    STEP    2 
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Figure  IV-3-3.   CASE  3,  STEP  3 


108 


*»««•*#*****«*«*•**•*«•«*«**«*£)««#»#**••»»««***«««»«**«*««*«# 


« 
« 

ft 

* 

ft 

* 

« 
« 
* 

♦ 
* 
* 
* 
* 
» 
* 
* 
* 

# 
* 
* 
* 
* 


0    0  0 
0  o 

0  0       0   00 
0    0    0  0 
0   0  00  0  00 
0   0    0 
0        8    0 
0  0   0      00 
0   00 
0   00   0 
0 
0 
0 


A 


0 
0  0 
0 
00 
0     0   0 
0      0 
080 
0 
0 
8       0 
0 
— 0 


« 
* 

* 
« 
* 
# 
* 

* 

« 

« 
# 


® 

X   X 


X      X  XX 
X   X 
XXX 


X 
X 


XX 
X    X 
XXX 
E  XX  E    X 
X 
XX   EX  X 
EXE     X 
X   E  X  X 
X    X     XX 

X    XX 
X  X 


« 
* 
ft 
♦ 
* 
* 
* 
* 

* 
# 
# 

ft 
# 
# 
* 


«###**■»#*#♦♦*#****#♦#*♦*♦*#  #**g#*-*#  ##*#♦♦#*#♦##»*  ####***#**»* 


Figure    IV-3-4.      CASE    3,    STEP    4 
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Figure    IV-3-5.      CASE    3,    STEP    5 
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Figure    IV-3-6P.      CASE    3,    STEP    7 
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Figure  IV-4-1.   CASE  4,  RUN  1,  STEP  1 
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Figure  IV-4-2.   CASE  4,  RUN  1,   STEP  2 
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Figure  IV-4-3.   CASE  4,  RUN  1,  STEP  4 
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Figure  IV-4-4.   CASE  4,  RUN  1,  STEP  5 
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Figure  IV-4-4P.   CASE  4,  RUN  1,  STEP  5 
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Figure    IV-4-5.      CASE    4,    RUN    1,    STEP    6 
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Figure  IV-4-6.   CASE  4,  RUN  2,  STEP  1 
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Figure  IV-4-7.   CASE  4,  RUN  2,  STEP  2 

A(l)  =  .1822    A(5)  =  .4875     B(l)  =  .6226  B(5)  = 

A(2)  =  .4815     A(6)  =  .2388    B(2)  =-.3740  B(6)  = 

A(3)  =  .3773    A(7)  =  .3191     B(3)  =-.0270  B(7)  = 

A(4)  =  .2770    A(8)  =  .3455     B(4)  =  .3068  B(8)  = 


-.3942 
.4340 
.1666 
.0788 
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Figure  IV-4-8.   CASE  4,  RUN  2,  STEP  3 

A(l)  =-.2417     A(5)  =  .6368     B(l)  =  .4122     B(5)  =  .1287 

A(2)  =-.0358     A(6)  *  .3585     B(l)  =  .6317     B(6)  =  .0277 

A(3)  =-.0299     A(7)  =  .4150     B(3)  =  .4963     B(7)  =  .0856 

A(4)  =-.0608     A(8)  =  .4790     B(4)  =  .3939     B(8)  =  .0702 
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Figure  IV-4-9.   CASE  4,  RUN  2,  STEP  4 

A(l)  =  .4496     A(5)  =  .0091     B(l)  =-.1616     B(5)  =  .6497 

A(2)  =  .6275     A(6)  =-.0388     B(2)  =  .0812     B(6)  =  .3574 

A(3)'=  .4934     A(7)  =  .0077     B(3)  =  .0621     B(7)  =  .4237 

A(4)  =  .3984     A(8)  =-.0193     B(4)  =  .0128     B(8)  =  .4837 
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Figure  IV-4-10.  CASE  4,  RUN  2,  STEP  6 

A(l)  =  .2460     A(5)  =-.6436  B(l)  =  .0580     B(5)  =  .1556 

A(2)  =-.2372     A(6)  =-.1104  B(2)  =-.8950     B(6)  =  .2618 

A(3)  =  .1528     A(7)  =-.3966  B(3)  =  .2344     B(7)  =  .1001 

A(4)  =  .2460     A(8)  =-.4644  B(4)  =-.1534     B(8)  =  .1202 
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Figure    IV-4-10P.      CASE    4,    RUN    2,    STEP    6 
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Figure    IV-5-1  .      CASE    5,    STEP    1 
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Figure  IV-5-2 
A(l)  =-.7748 
A(2)  =  .1863 
A(3)  =-.1346 
A(4)  =  .5890 


CASE  5,  STEP  2 
B(l)  =-.2831 
B(2)  =-.6366 
B(3)  =  .7173 
B(4)  =-.0072 
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Figure  IV-5-3. 
A(l)  =  .1429 
A(2)  =  .5127 
A(3)  =-.8460 
A(4)  =  .0318 


CASE  5,  STEP  3 

B(l)  =-.8365 

B(2)  =-.0491 

B(3)  =-.1513 

B(4)  =  .5244 
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Figure  I V- 5-4 
A(l)  =-.0571 
A(2)  =  .1522 
A(3)  =-.8827 
A(4)  =-.4409 


CASE  5,  STEP  5 

B(l)  =  .8120 

B(2)  =  .5511 

B(3)  =-.0502 

B(4)  =  .1856 
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Figure    IV-5-4P,    CASE    5,    STEP  5. 
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Figure  IV- 5-5. 
A(l)  =  .4220 
A(2)  =  .3291 
A(3)  =-.7281 
A(4)  =-.4284 


CASE  5,  STEP  6 

B(l)  =-.8145 

B(2)  =-.2528 

B(3)  =-.5022 

B(4)  =-.1429 
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Figure  IV-5-6. 
A(l)  =  .0767 
A(2)  =  .3515 
A(3)  =-.8091 
A(4)  =-.4647 


CASE  5,  STEP  12 

B(l)  =  .9772 

B(2)  =  .1482 

B(3)  =  .1515 

B(4)  =  .0095 
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Figure  IV-6-6.  CASE  6,  RUN  2,  STEP  12 
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Figure  IV-6-7.   CASE  6,  RUN  2,  STEP  13 
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Figure  IV-7-3.   CASE  7,  RUN  1,  STEP  11 
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Figure  IV-7-5.   CASE  7,  RUN  1,  STEP  20 
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Figure  IV-7-8.   CASE  7,  RUN  2,  STEP  2 
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Figure  IV-7-9.   CASE  7,  RUN  2,  STEP  9 
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Figure  IV-7-10.   CASE  7,  RUN  2,  STEP  12 
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Figure  IV-7-11.   CASE  7,  RUN  2,  STEP  14 
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Figure  IV-7-12.   CASE  7,  RUN  2,  STEP  27 
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